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For equal-spin s\ tetramer single molecule magnets with ionic point groups g = Tj, D^, D2d, 
Cih, Ci v , and Sa, we write the group-invariant single-ion, Dzaloshinskii-Moriya, and anisotropic 
near-neighbor and next-nearest-neighbor exchange Hamiltonians using the respective local axial and 
azimuthal vector groups. In the molecular representation, the symmetric anisotropic exchange in- 
teractions renormalize the near-neighbor and next-nearest-neighbor isotropic exchange interactions 
to J g and J' g . The local single-ion and symmetric anisotropic exchange interactions generate the 
site-independent anisotropy interactions in the molecular representation, Jf (/if ) , Jf z (/xf 2 ) , and Jf z , 
respectively, where \x\ and represent the group-consistent sets of parameters required to diago- 
nalize the single-ion and symmetric near-neighbor and next-nearest-neighbor anisotropic exchange 
Hamiltonians, respectively. The site-independent group-invariant near-neighbor and next-nearest- 
neighbor Dzaloshinskii-Moriya interactions generally may be written as vectors of rank three and 
two, d 9 and d'j 3 , respectively. The local anisotropy interactions also generate site-dependent single- 
ion, near-neighbor, and next-nearest-neighbor anisotropic exchange interactions in the molecular 
representation. Using our exact, compact forms for the single-ion matrix elements, we evaluate the 
eigenstate energies to first order in Jf (juf), Jf z (/^ 2 )j aQ d Jz,z- There are two types of ferromagnetic 
(FM) ( J g > and antiferromagnetic (AFM) (J g < tetramers. In Type I, J' g — J g > 0, the tetramer 
acts as two dimers with the maximal pair quantum numbers S13 = S24 = 2si at low temperatures 
T. Type II tetramers with J' g — J g < are frustrated, with minimal values of the pair quantum 
numbers S13 and S24 at low T. For both Type I and II AFM tetramers, we evaluate the first-order 
level-crossing inductions analytically for arbitrary si, and illustrate the results for si = 1/2, 1,3/2. 
Accurate Hartree expressions for the thermodynamics, electron paramagnetic resonance (EPR) ab- 
sorption and inelastic neutron scattering cross-section are given. A procedure to extract the effective 
microscopic parameters for Types I and II FM tetramers using EPR is given. 

PACS numbers: 75.75.+a, 75.50.Xx, 73.22.Lp, 75.30.Gw, 75.10.Jm 



I. INTRODUCTION 

Single molecule magnets (SMM's) have been a topic of 
great interest for more than a decade, because of their 
potential uses in quantum computing and/or magnetic 
storage, which are possible due to magnetic quantum 
tunneling (MQT) and entangled states. In fits to a wealth 
of data, the Hamiltonian within an SMM was assumed 
to be the Heisenberg exchange interaction plus weaker 
total (global, or giant) spin anisotropy interactions, with 
a fixed overall total spin quantum number s. MQT and 
entanglement were only studied in this simple model. 

The simplest SMM's are dimers. 0,0] Surprisingly, two 
antiferromagnetic dimers, an Fe2 and a N12, appear to 
have substantial single-ion anisotropy without any appre- 
ciable total spin anisotropy. |4], |fj Although the most 
common SMM's have ferromagnetic (FM) intramolecu- 
lar interactions and contain n > 8 magnetic ions,0 a 
number of intermediate-sized FM SMM's with n = 4 and 
rather simple molecular structures were recently studied. 
The Cu 4 tetramer Cu 4 OCl 6 (TPPO) 4 , where TPPO is 
triphenylphosphine oxide, has four spin 1/2 ions on 
the corners of a regular tetrahedron, with an s = 2 
ground state and approximate Td symmetry. [ToL ITll Il2^ 



In this case, there are no single-ion anisotropy effects, 
but anisotropic symmetric exchange interactions were 
thought to be responsible for the zero-field energy 
splittings. [H The Co 4 , Co4(hmp) 4 (MeOH) 4 Cl4, 
where hmp is hydroxymethylpyridyl, and Cr 4 , 
[Cr 4 S(0 2 CCH 3 ) 8 (H 2 0) 4 ](N03)2-H20, compounds 
have s = 6 ground states with spin 3/2 ions on the 
corners of tetrahedrons. |{j 1 1 1| Those compounds have S4 
and approximate Z?2d symmetry, respectively. 0, ^| A 
number of high symmetry s = 4 ground state Ni 4 struc- 
tures with spin 1 ions were reported. [IJI, [lj| [lT], lla liil 
Two of these, [Ni(hmp)(ROH)Cl] 4 , where R is an alkyl 
group, such as methyl, ethyl, or 3,3-dimethyl-l-butyl 

form tetramers 
Two others, 

Ni 4 (ROH)L 4 , where R is methyl or ethyl and H2L 
is salicylidene-2-ethanolamine, had approximate S4 
symmetry, although the precise symmetry was only 
Ci.[l5( Several planar Mn 4 compounds with the Mn +3 
spin 2 ions on the corners of squares were made, with 
overall s — 8 tetramer ground states. |io| Although two 
of these complexes had only approximate S4 symmetry, 
one of these complexes, Mn 4 Cl 4 (L') 4 , where H2L' 
is 4-i-butyl-salicylidene-2-ethanolamine, had perfect 
S4 symmetry. [20j Inelastic neutron scattering (INS) 



and hmp is 2- hydroxymethylpyridyl 
with precise 5 4 group symmetry. |l8l |l 
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FIG. 1: T d (left) and D 4h (right) ion sites (filled). Circle: 
origin. Arrows: local axial z^ d (left), azimuthal x^ 4h (right) 
single-ion vectors. The axial vectors z® ih = z, normal to the 
ionic plane. 



FIG. 2: D 2d (left) and S 4 (right) ion sites (filled). Circle: ori- 
gin. Arrows: local axial single-ion vectors. The g — D2d,S 4 
axial vectors zf make the angles 



ft axial vector if 4 also makes the angle a\ with the x axis, 
where cosqi = sin#f 4 cos</>f 4 . 



f with the z axis, and the 



experiments provided strong evidence for single-ion 
anisotropy in C04 and a N14 with approximate S 4 
symmetry. 0, The presence of single- ion or exchange 
anisotropy actually precludes the total spin s from 
begin a good quantum number. Q Fits to electron 
paramagnetic resonance (EPR) Ni 4 data assuming a 
fixed s were also problematic, suggesting single- ion or 
exchange anisotropy in that tetramer, as well.|2l| 

Recently there have been microscopic treatments of 
dimers,Q trimers, and tetramers, including Zeeman g- 
tensor anisotropy, single-ion anisotropy, and anisotropic 
exchange interactions. |22j Most of those treatments and 
their recent extensions to more general systems expressed 
the single-spin matrix elements only in terms of Wigncr 
3j, 6j, and 9j symbols. [22I |23| While such treatments 
are very helpful in fitting experimental data, more com- 
pact analytic forms are desirable to study microscopic 
models of FM SMM's in which the MQT and entangle- 
ment issues crucial for quantum computing can be un- 
derstood. We constructed the single-ion and anisotropic 
near-neighbor (NN) and next-nearest-neighbor (NNN) 
exchange SMM Hamiltonians from the respective local 
axial and azimuthal vector groups for equal-spin tetramer 
SMM's with point group symmetries g = Td, D 4 h, Did, 
Cih, C 4v , and S 4 , and found compact analytic expres- 
sions for the single-spin matrix elements of four gen- 
eral spins. Surprisingly, each local vector group gen- 
erates site-dependent molecular single-ion and exchange 
anisotropy. We evaluate the magnetization, specific heat, 
EPR and INS transitions in the Hartree approximation, 
and provide a procedure for extracting the effective site- 
independent microscopic parameters using EPR. 



II. STRUCTURES AND BARE HAMILTONIAN 

For SMM's with ionic site point groups g = Td, D 4 h, we 
assume the four equal-spin s\ ions sit on opposite corners 
of a cube or square of side a centered at the origin, as 





FIG. 3: C 4h (left) and C 4v (right) ion sites (filled). Circle: 
origin. Arrows: local azimuthal x^ ih (left) axial z% 4qj (right) 
single-ion vectors. The axial vectors z% 4h — z. The C 4v axial 
vectors each make the angle 81 = tv/2 — 8[ with with z axis. 
The dotted arrows (equivalent to the D 4 h azimuthal single-ion 
are their projections in the xy plane. 



vectors 



pictured in Fig. 1. For SMM's with g — D2d,Si, we 
take the ions to sit on opposite corners of a tetragonal 
prism with sides (a, a, c) centered at the origin, as in Fig. 
2. The ions for g — 0^,0^ also sit on the corners of 
a square of side a centered at the origin, as pictured in 
Fig. 3, but the ligand groups have different symmetries 
than for the simpler D 4 h case pictured in Fig. 1 . In each 
case, we take the origin to be at the geometric center, so 
that Y^ n =i Vn = 0' where the relative ion site vectors 



r n = |[-7^ + 7nW]-|(-l) B *. 
7± = 4(-l)" /2 ±^(-l) (n+1)/2 , 



(1) 
(2) 



and e± = [1 ± (-1)"]/2.[2J For clarity, the values of 
and 7^ are given in Table 1. In tetrahedrons with 
g = Td, c/a = 1, approximately as in Cu4.[l0| In squares 
with g = D 4 h, C 4 h, or C 4v , c = 0. The high D 4 h symme- 
try is approximately exhibited by the square Nd 4 com- 
pound, Nd4(OR)i2, where R is 2,2-dimethyl-l-propyl, in 
which the Nd +3 ions have equal total angular momen- 
tum j = 9/2. [HHl] We note that the Mn 4 SMM's with 
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^ r fn 'Yn 

10 1-11 

2 1 0-1-1 

3 1 1-1 

4 10 11 

TABLE I: Values of e± and 7* 



approximate or exact S 4 symmetry also have c = O.[20j 
In tetragonal prisms with g = D 2 d or S4, c/a > 1, ap- 
proximately as in Co4,|b| or c/a < 1, as in Mn4 and a 
Nu . fl6t Eof x, y, z are the molecular (or laboratory) axes 
of each SMM. 

The most general Hamiltonian quadratic in the four 
spin operators S n may be written for group g as 



H 9 



J2B-g g n -S n + J2 Sn-r> a n>n ,-s n ,, (3) 



where [lb is the Bohr magneton and B = 
B (sm 8 cos cf>, sin 8 sin cf>, cos 8) is the magnetic induc- 
tion at an arbitrary direction (8, <£>) relative to the 

molecular coordinates (x. v. £1.|22Ll27j 

The exchange matrices D 9 , for n! 7^ n can have an- 
tisymmetric terms, which are the Dzyaloshinshii-Moriya 
(DM) interactions. Of the six molecular group symme- 
tries under study, the DM interactions are non- vanishing 
for all but the highest group symmetry, 7d.|28|,|29j For the 
other five group symmetries, the DM interactions lead to 
weak perturbations to the eigenstate energies. Inequiva- 
lent ions or ligands, or local structural distortions from 
the ideal point group symmetry, as commonly occur in 
real systems, |30j break these g symmetries, introducing 
stronger DM interactions. Such lower symmetry groups 
will be discussed elsewhere. |3ll) 

For simplicity, we take *g 9 n to be diagonal, isotropic, 
and site-independent, so that the Zeeman interaction 
may be written in terms of a single gyromagnetic ra- 
tio 7 ~ 2/ib- Thus in the following, g only refers to the 
molecular group. We separate D 9 n , into its symmetric 
and antisymmetric parts, D 9 nn , — D 9 ^ s n , + D 9 ^ a n ,, re- 
spectively. For n' — n, the single-ion D 9 L n is necessarily 

symmetric, so D^ a n = 0. For each g, the four D^' s n con- 
tain the local single-ion structural information, and the 
six distinct symmetric D„ n i contain the local symmet- 
ric exchange structural information, which lead to the 
isotropic, or Heisenberg, exchange interactions, and the 
remaining symmetric anisotropic exchange interactions. 
The six distinct antisymmetric D 3 ' , contain additional 
local structural information which lead to the DM inter- 
actions. Physically, the symmetric anisotropic exchange 
interactions also contain the dipole-dipole interactions, 
which can be even larger in magnitude than the terms 
originating from actual anisotropic exchange. |27j 



As is well known, each of the symmetric rank-three 
tensors (or matrices) can be diagonalized by three rota- 
tions: a rotation by the angle (j) 9 n , about the molecular z 
axis, then a rotation by the angle 8 9 n , about the rotated 
x axis, followed by a rotation by the angle ipn a bout 
the rotated z axis, This necessarily leads to the three 
principal axes x n n ,, y n n ,, and z n n ,. For the single-ion 
axes with n' = n, we denote these principal axes to be 
x n , y n , and z n , respectively, which are written explic- 
itly in Sec. III. The non-vanishing matrix elements in 
these locally-diagonalized symmetric matrix coordinates 

are ^n,n',»' ^n,n', yy and J^n',..- SinCe the Structural 

information in each of the D 9 ^ s n , depends upon the local 
environment, in the absence of molecular group g symme- 
try, each of these angles would in principle be different 
from one another. We characterize the single-ion and 
symmetric anisotropic exchange vector sets by those for 
the spin at site 1, the NN pair at sites 1 and 2, and the 
NNN pair at sites 1 and 3, respectively, 



W 3 > 



(4) 
(5) 
(6) 



respectively, where 8f = 8\ l , 8\ 2 = 8\ 2 , and 8\ 3 = 8\ 3 . 

Since the antisymmetric exchange matrix cannot gen- 
erally be diagonalized, instead of writing a different gen- 
eral set of exchange vectors, we first employ the Moriya 
rules for the interactions between each pair of spins in 
the molecular representation, 27, 29] and then impose the 
group operations in that representation. This leads to an 
additional set of as many as five parameters for the NN 
and NNN DM interactions, which we write as the three- 
and two- vectors, d 9 and d_[ 9 , respectively. 

For the six high-symmetry groups under study, we ana- 
lyze the effects of molecular group symmetry upon these 
coordinates. As shown in the following, the molecular 
group symmetry can greatly simplify the theoretical anal- 
ysis of the measurable quantities, so that only the var- 
ious interaction magnitudes and these five angular sets 
are required. Some group symmetries impose severe re- 
strictions upon the elements of the angular sets. 

In the absence of any anisotropy interactions, the bare 
Hamiltonian 7Yq is given by the Zeeman and Heisenberg 
interactions, 

n 9 = - 7 b • s - j' g {Si -s 3 + s 2 - s 4 ) 

-J g (S 1 -S 2 + S 2 -S 3 + S 3 -S 4 + S 4 - Sx), (7) 
which can be rewritten as 

where we have dropped the irrelevant constants, and 

J' Td = M, (9) 
J ' 9 * J 'y ( 10 ) 
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for g = D 2 d, S4, D 4 h, C4/1, and C 4v . In terms of the diag- 
onalized matrix elements, — 2J„ 



D g,s 



D x%,m and 



-2J' g = Dff 3xx + Dfll yy , for instance. 



We note that for c/a < 1, — J g and — J' g are the NN and 
NNN Heisenberg interactions for D 2 d, S4, as for Ci V , C 4 h 
and D^h with c = 0, but for c/a > 1, — J g and — are 
the NNN and NN Heisenberg interactions, respectively, 
for D 2d and S4. In Eq. JHJ, S13 = S1 + S3, S 2 4 = S2+S4, 
and S — S13 + S24 is the total spin operator. 26] 



III. THE SINGLE-ION AND ANISOTROPIC 
EXCHANGE HAMILTONIANS 



In order to properly take account of the molecular 
group g symmetries, it is useful to describe the single- 
ion and symmetric anisotropic exchange interactions in 
terms of the local coordinates. We write the local Hamil- 
tonian for these interactions, and require that it must be 
invariant under all A operations 9 x of g. As noted above, 
we write the antisymmetric exchange interactions in the 
molecular representation, employing the group symme- 
try interactions relevant to the particular pair involved 
in each exchange, and then impose the overall group sym- 
metry. In Sees. IV-IX, we then impose the group sym- 
metries on these interactions for C^h, D4h,C'4 V , S4, D 2 d, 
and Td molecular group symmetries, respectively. 



A. Local single-ion Hamiltonian 

To do so for the single-ion anisotropy, we first define 
the vector bases for each group g. We define the local 
axial vector basis to be {z n } for each g, and the orthog- 
onal azimuthal vector bases to be {x^} and {y n }- The 
elements of these bases are the vectors that diagonalize 

the single ion matrix from D 9 n n to D 9 n n . Since we 
employ these vectors repeatedly, we write them here for 
simplicity of presentation. The diagonalized vectors x n , 
etc., may be written in the molecular (x,y,z) represen- 
tation as 



X ... 



~9 
Vn 



a 9 



cos <§fl n cosi/4 



cos ip 9 ^ sin cj) 9 n 



n 

)■" 

it 

sin d 9 n 



cos 9 a n 
. cos. 7 , 
9 n sin^ 



sin (p 9 ^ sin ^ 
9 9 n cos cj> 9 n sin ip\ 



- cos (p 9 n sin tp^ — cos Q 9 n si 

- sin ip n sin (p 9 n + cos d 9 n c< 

sin Q 9 n cos 



cos 0% sin (p 9 n cos tp n 
cos <p 9 n cos ip 9 



b 9 

rn 



sin0£sinv„ 
- sin Q 9 n cos 4> 9 n 
cos 6 9 




which satisfy x n x y n = z n . We then write the most 
general quadratic single-ion anisotropy interaction as 



n 9 

1 i- o 



^ 1 (ja,n(Sn ' % n ) 



+J 9 J(S n -x 9 n f-(S n -y:f}), (14) 

in terms of the site-dependent axial and azimuthal in- 
teractions Ja n , analogous in notation to that for 
homoionic dimers.jj] In terms of the diagonalized matrix 
elements, - J| >n = D 9 n ^ zz - (D^ xx + D 9 n % m )/2 and 
— J 9 =(f) 9 < s — f) 9 < s )/2 

We require H, 9 s l to be invariant under all allowed g 
symmetries. As we shall see for the six high-symmetry 
cases under study, this forces J 9 n and J|„ to be in- 
dependent of n for each g in these local coordinates, 
and places constraints upon the single-ion parameter set 

{0 9 n,KM- 



~. 9 \2- 



Local symmetric anisotropic exchange 
Hamiltonian 



In addition to the single-ion interactions, the other 
microscopic anisotropic interactions are the anisotropic 
exchange interactions, which include the dipole-dipole 
interactions. We neglect intermolecular dipole-dipole 
interactions, which can be important at very low 
temperatures, [s^j As for the single-ion interactions, 
we first construct the symmetric anisotropic exchange 
Hamiltonian TL 9 ae in the local group coordinates. In this 
case, there are distinct local vector sets for the NN and 
NNN exchange interactions. Diagonalization of the sym- 
metric anisotropic exchange matrix D^ a n i leads to D 9 n s n , 

and the vector basis {x nn t,y nn ,,z nn ,}, given by Eqs. 
(|1 1(1 - 1(13(1 with the subscript n replaced by n,n'. 

The local symmetric anisotropic exchange Hamiltonian 
H%£ is then generally given by 



2 6— 2m 



U 9,l 



'y ] ^ ] \jn,n+m{^n ' n n+m ){S n + m ■ Z n n+m ) 



m=l n=l 



~(S n ■ y n n+m )(S n + m ■ y n n + m ) 



(15) 



where we define S$ = Si, as if the four NN spins were on 
a ring. In Eq. 1)15(1 . the axial and azimuthal interaction 
strengths = D 9 n %,. zz -{D 9 n ^. xx + D 9 n %,. yy )/2 and 

- J ni> = " 5 n,n', w )/ 2 ' as for the single-ion 

interaction strengths, but in this case n' — n + m, where 
to = 1,2. 

Since the structures allow for two distinct bond 
lengths, with four NN's and 2 NNN's for c/a < 1, there 
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are two sets of local axial and azimuthal coordinates. 
The vector sets {I^ >n+1 }, {x^ n+1 }, and {y 9 n<n+1 } each 
have four elements, and are the bases for the NN ax- 
ial and azimuthal symmetric anisotropic exchange vec- 
tors. In addition, the vector sets {z nn+2 }, {x nn+2 }, 
and {if n n+2 } each have two elements, and are the NNN 
symmetric anisotropic exchange vector bases. These are 
both different from the single-ion vector bases. 

For the six high-symmetry groups under study, the 
groups symmetries force the interaction strengths Jnn+m 
and Jnn+m to be independent of n, but are generally dif- 
ferent for NN (to — 1) and NNN (m = 2) interactions. 
The group symmetries also place restrictions upon the 
NN and NNN symmetric anisotropic exchange param- 
eter SetS {0 9 n t n+m^ 9 n+ra^n,n+m} for m = 2 - These 

restrictions are much stronger for the NNN parameter 
set than for the NN parameter set. 



place restrictions upon the d 9 n n+m , causing them to be 
independent of n, so that only the two three-vectors d\ 2 



'13 



tetramer group. However, some components of the DM 
interaction vectors can have site-dependent signs. 

In Sections IV-IX, we evaluate H 9 si , Ti. 9 e , and H 9 DM for 
g = C ih ,D ih , C iv ,S4, D 2d , and T d , respectively. 



IV. C ih GROUP SYMMETRY 

We first discuss the simplest group symmetries, C^h- 
Besides the trivial identity operation, the allowed group 
operations O x 4h for A = 1,2,3 are clockwise and coun- 
terclockwise rotations by n/2 about the z axis, and re- 
flections in the xy plane. [24j These operations are repre- 
sented respectively by the matrices 



C. Antisymmetric anisotropic exchange Hamiltonian 

As noted above, since the antisymmetric anisotropic 
exchange matrix with three real distinct matrix elements 
cannot generally be diagonalized by a unitary transfor- 
mation, there is no particular vector basis resulting from 
the diagonalization. We therefore write the antisymmet- 
ric anisotropic exchange, or DM, Hamiltonian Ti. 9 DM in 
the molecular representation. |27l| 



V 9 

rL DM 



E 



6-2m 

E 

n=l 



d 9 

n,n-\-m 



(s n x S n+m \ . (16) 



We note that in these molecular coordinates, the DM 
interaction three-vectors d 9 n n + m depend explicitly upon 
the exchange bond indices n,n + m for each group g. We 
then employ the local group symmetries to relate them 
to one another. 

The rules for the directions of the d 9 n n + m were given 
by Moriva.|29| and were employed for a dimer example 
by Bencini and Gatteschi,|27j The Moriya rules are: (1) 
d 9 n , vanishes if a center of inversion connects r n and 
tv. (2) When a mirror plane contains r„ and tv, d 9 n n , 
is normal to the mirror plane. (3) When a mirror plane 
is the perpendicular bisector of r„ — r n > , d 9 n n , lies in the 
mirror plane. (4) When a two-fold rotation axis is the 
perpendicular bisector of r n — tv , then d 9 n , is orthogo- 
nal to r n — tv • (5) When r n — TV is a p-fold rotation axis 
with p > 2, then d 9 n n , is parallel to r n — TV- As noted 
above, we shall incorporate these rules in the molecular 
representation. 

In the following, we impose the Moriya rules on each 
anisotropic exchange pair, and then impose the required 
group symmetries on the six pairs. For each g, the sym- 
metries relate the NN and NNN DM interactions to one 
another, greatly restricting the number of parameters. 
For the six groups under study, the group symmetries 



1,2 



Of 4 " = 




(17) 



(18) 



single-ion anisotropy 



We require 0? Ah H% 



= U^ h for A = 1, 2, 3. We first 
consider the effects of Of 4h . We note that Of 4h r„ = 
r n +i- Hence, this operation rotates the crystal by 7r/4 
about the z axis. Graphically, this is accomplished by 
relabelling Fig. 3(a) with n — ► n — 1, so that S\ is now 
at r 2 , etc. Thus, 



0^ h H 



4 

E 

n=l 



( J C ^(S . z Clh \ 2 

\ J a,n \ J n ^n+l) 



+ ^ejn' 1 [{&n ' ^n+l) {&n ' Vn+l) . 



(19) 



z n+ i, etc. In order that the axial 
part of H. si is invariant under the transformation, we 
require 



where 0^ ih z^ h 



Ja 



and 



sin 0n Ah sin <^ 4 " 
sin6>£ 4h cos0^ 4 ' 1 
cos 6% ih 



-sin^cos^!, 



(20) 



(21) 
(22) 
(23) 
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There are two possible solutions for these recursion rela- 
tions. One possibility is 



9 C4h = 7T 

b^! 4h = arbitrary. 



The second possible solution is 

u n — u n-l> 



(24) 
(25) 



(26) 
(27) 



However, imposition of the xy mirror plane symmetry, 



0^ 4h , such as setting 0^ 4h z^ 4h = ±z^ 4h forces either 
the even parity 9% 4h = or the odd parity 9^ 4h = 1 . We 
choose the odd parity solution, setting 



6? 4h = 0, 



(28) 



which implies z n = z, the axis of high symmetry. 

Now, we also must also require the azimuthal part of 
H° 4h ' e to be invariant under Of 4fe . Using Eq. ©tft, we 



see from Eqs. and (|12() that x n — ysmx 

xcosXn 4h i etc, where 



X«*" = ^ 4h +^ 4h - (29) 

en c 
Then, invariance of under 0\ 4h requires 

Je,ri = Je,n-1 = ^& (30) 

and 



cosx£ 4h = sin X ^\, 



sin x 



C4/1 



that 



An An— 1 



(31) 
(32) 

(33) 



We note that for these vector groups, T~t si 4h ' is automat- 
ically invariant under 2 4 and C 3 4 . In the latter case, 

Of 4 i^ 4h = — £^ 4h has odd parity. 

Since Of 4h interchanges the x and y components of 
the azimuthal vectors, it is useful to define the functions 

<W) = e+(-l)" /2 sin0 + e-(-l)(" +1 )/ 2 cos0,(34) 
<TO) = 4(-l)" /2 cos^-e;(-l)(" +1 )/ 2 sin0.(35) 

We note that <5±(tt/4) = 7^/%/2 and S±(3tt/4) = 
±7^/v2- For clarity, these new functions are listed in 
Table II. Some useful relations are 

(C(^)) 2 +( ( 5-(0)) 2 = 1, (36) 

(CW) 2 -(W)) 2 = (-l)"co S (20), (37) 

= ^7^sin(2^. (38) 



n 8t{4>) 5 n (<ft) 

1 — cos </> sin (f> 

2 — sin — cos 

3 cos 4> — sin 

4 sin <f> cos 

TABLE II: Values of the functions 5± (</>). 



Hence, the single-ion vector groups for g = C^h sym- 
metry are {x^ 4h }, {Vn ih }i {^n 4h }' where the elements 
are given by 



-e(xi 4h )*+<uxr i )y, (39) 

-5-(x? 4h )x-S+( x ? 4h )y, (40) 



z = z. 



The single-ion parameter set for g = C^h is 

x ° 4h = 0f 4h + vf". 



(41) 



(42) 
(43) 



C4h symmetric anisotropic exchange 

The construction of the symmetric anisotropic ex- 
change vector groups is entirely analogous to those of 
the single-ion anisotropy. We find 



rf,Cih _ jn— l,n+m— 1 _ j 

J n,n+m ~ J f ~ J f,m> 

— Jfj'tn? 

jc,Cih Tn— l,n+m — 1 

^n,n+m ^ c 

= Jc.mi 



(44) 
(45) 
(46) 
(47) 



and the near-neighbor and next-nearest-neighbor vector 
groups are 



aC+i = Si(Xii h )x + 6-(xi 2 ih )y, (48) 

yS+l = - 6 n(X?2 ih )x-tn(X?2 ih )y, (49) 



z 



and 



n,n+l 



*^n,n+2 
2/n,n+2 



(50) 



-S-(x? 3 ih )x-S+( x ?i h )y, (52) 



(53) 



respectively, where 



(54) 



for p — 2,3. 
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C. dh antisymmetric exchange interactions 

With antisymmetric exchange, we first impose the 
Moriya rules for C 4 h symmetry. We first note that 



ox 4h o\ 



-o 



3 ) 



of inversion for the next-nearest neighbor bonds with 
g = CwHi This implies 



(55) 



Since all four spins lie in the xy mirror plane, C 3 4h , 
we have 



d Cih — rf 
"n,n+l — 



(56) 



Then, imposing the O l 4h rotational symmetry on H D 4 ^ , 
we find 



"n.n+l — u n-l,n 

= d z 



(57) 
(58) 



Hence, we require (D? ih xia — ±X2 4, 0? 4h xi^ = ±£43, 



0° ih xi, 3 = ±xi, 3 , 0^ h x 2A = ±x 4 ,2, and 0^ fc x 1)3 = 
±£3,1, Of 4, *&2,4 = ±^2,4- The even (odd) parity choices 
lead to sin \\ ih = icosxi 4 " 1 , and we make the even 



so that the origin is a center parity choice, 



r> 4 
Xi 



7T/4, 



(62) 



consistent with Fig. 1. We then obtain the single-ion 
vector sets for D 4 h symmetry, 



1 1 +- 

-7=2^ 

z, 

(0.7T/4), 



InV) 



(63) 

(64) 

(65) 
(66) 



where the second entry refers to X\ ih = ( t ) i ih + 4>i ih 
7r/4. We note that x^ 4h — r n , as pictured in Fig. 1. 



is a scalar independent of n. Thus, the antisymmet- 
ric anisotropic exchange interaction for C 4 h symmetry 
is simply 



K 



DM 



d z ^z ■ (S n x S n+ i) 



(59) 



V. D 4h GROUP SYMMETRY 

We now consider the higher symmetry g — D 4 h- This 
group contains the three symmetries of C 4 h, so that 

QD ih = a C ih for A = 3 In addition; thcre are four 

two-fold rotations about axes normal to the principle z 
axis. These rotation axes are the x and y axes and the 
y = ±x diagonals. These D 4 h operations are represented 
by 



4,5 







( ±1 

T-l 
y 0-1 

±1 



±1 






-1 



(60) 



(61) 



A. D 4 h single-ion anisotropy 



We begin with the same invariances of Ti^ 4 ' 1 ''' as for 
Cih symmetry, leading to the single-ion vector groups 
given by Eqs. with C 4 h —> D 4h . As for C 4 h 



symmetry, = J°%_ x = J a and = J D 



e,n— 1 



D 4 



J e . We then impose the additional symmetries O 
for A = 4 — 7. We note that OP 4h ri = n, for example. 



B. Da,h symmetric anisotropic exchange 

As for the single-ion interactions, we begin with the 
symmetric anisotropic exchange vector groups for C 4 h- 



As for Cih symmetry, we have Jn\n+m = J\ 

J c ,m ana J n \ n+m — •J n -l,n+m~l 



c,D ih 
'n— l,n+m— 1 



Jf,m ■ We then impose 



J 1\ 



the additional symmetries, setting 4 x 12 

— £C 12 and 4 x 13 h = x 24 , for examples. Thus, it is 
elementary to obtain the restrictions 



cosxf 2 4i 
cosxfs" 



sin X f 2 4 ' 



and we choose the solutions 

A12 



A13 — 



7T 

2' 

7T 



(67) 
(68) 

(69) 
(70) 



These solutions are consistent with the remaining sym- 
metries O® 4 for A = 5, 6, 7. Thus, we obtain the 
near-neighbor and next-nearest-neighbor D 4 h symmetric 
anisotropic exchange vector groups, 



X n,n+1 



Vn,n+X 
Un,n+1 



1 



1 



(in + ln) X + o(7n 



Alt ~ln) x ~ o(7n + In 



and 



"n,n+2 



Un,n+2 



Z n,n+2 



V2 



(-lt X + lnil), 



V2 



(ln X + lnV), 



(71) 
(72) 

(73) 

(74) 
(75) 
(76) 
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respectively. 



C. Dih antisymmetric anisotropic exchange 

For D^h symmetry, the origin is a center of inversion 
for the next-nearest-neighbor DM interactions, forcing 
them to vanish, as for C^h- Moreover, it is easy to show 
that tetramers with D^h symmetry, the near-neighbor ex- 
change given by Eq. I|59|) also in invariant under 0® ih . 
Hence, the most general antisymmetric anisotropic ex- 
change Hamiltonian for D^h symmetry is 



(77) 



precisely the same as for C^h symmetry. 



VI. C iv GROUP SYMMETRY 

For d v , the group operations 0^ iv for A = 1 , . . . , 6 are 
clockwise and counterclockwise rotations by 7r/2 about 
the z axis (equivalent to G 1 f in Eq. ^T7\ ). and reflec- 
tions in the xz, yz, and diagonal mirror planes containing 
the z axis and the lines y = ±a;.[24| The additional op- 
erations are represented respectively by the matrices 



^3,4 



O 



5,6 




(78) 



(79) 



where the upper signs refer to the lower subscript num- 
bers. 



A. Ci v single-ion anisotropy 



which is satisfied by either choosing sin#^ 4 " = or 
sin cj>f Av = — cos <^ 4 " . The former solution is appro- 
priate for the higher-symmetry C'ih , D^h cases in which 
the xy plane is a molecular mirror plane, pictured in 
the right panel of Fig. 1, but this is not required for 
the most general tetramers exhibiting C^ v symmetry, for 
which the ligands on opposite sides of the ring may be 
different. We therefore choose the more general case, 
setting </)f 4 " = 37r/4, and leaving 6*f 4 " arbitrary. We 



-c 4v 



-c iv 



then construct z 3 *" and z i ™ from O x 4 " z 



-Civ 



--C 4 „ 



and 



Ci, 



-Ci v 



-Civ 



-Civ 



The resulting vector group {z n "} 



is invariant under the diagonal mirror plane operations 

c 4« „o n c *» i.° 4v — z C4v Q°iv ~ C4 " 



Ot%\ as 5 



O, 



c iv 



-Civ 



-Civ 



O, 



Civ 



-Civ 

Z l,3 



-1,3 



-Civ 

Z 3,l 



~ *1,3 > and 

Hence, these 



diagonal mirror planes introduce no restrictions upon 
6*p 4 ". We also note that the yz mirror plane operations 

„ %Civ ^Civ -Civ -Civ 

^i z i 4 — z i i ana <u^Z2 3 = z 3 2 ■ 



Similarly, from x 1 4 " given by Eq. with n = 1, g 



C\ v , we construct x 2 = 1 iv x 1 = 2 4v x 1 . After 



imposing our choice f^ 4 " = 37r/4, we find that these two 
equations are only satisfied if cos ipf 4 ' 
choose Vf 4 " = -tt/2, so that xf 4 " = 
as for D^h symmetry. 



= 0. We therefore 

Civ _1_ J,Civ 



+ =7r /4, 



By successive applications of Of 411 

-Ci 



and 3 iv , we ob- 



tain the three vector groups {x 
with the elements 



-Civ 



-Civ 

Vn 



-Civ 

z„ 




} and {aT"}, 



(81) 












V J 




-it sin 0f 4 


7-sin^ 4 " 


V^cosflF 4 " 



(82) 



(83) 



As for dh and Dih symmetries, we first impose the 
O^ 4 " invariance, the rotation by tt/2 about the z axis. 
This again leads to jf 4 " = J°n-i = J a and J e C4 " = 
J^n-i = Jet independent of n. As for the axial single- 
ion anisotropy with C^h symmetry, this leads to Eqs. 
(|2i p -l|23 p . with the two choices for the recursion rela- 
tions for #^ 4 " and <^„" . However, with G^ v symmetry, 
the xy plane is not a mirror plane. Instead, we im- 
pose the xz mirror plane operation O^ 4 " , and require 



±z x *" . Choosing the even parity case leads and ip^ 4v , leading to 



where the 7^ are given in Eq. © and Table I. These 
forms are also consistent with the yz and diagonal mirror 
plane operations (O x 4 " operations for A = 4,5,6), and 

leave Eq. Ijl4(l invariant. [24j We note that some of the 

-c^ 

reflections of some of the y n have odd parity. We also 
note that these vectors reduce to those for D±h symmetry, 
Eqs. dHH|- dnni> in the 6>f 4 " -> limit. We reiterate that 
the Ci v group symmetry imposes restrictions upon (j> 1 4 " 



to the equation 



sin £»f 4 " sin cf>f iv 



-sin(9f 4 " cos^f 4 ", (80) 



R 4 "} = 



9 Civ 



,37r/4,-7r/2). 



(84) 
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B. d v symmetric anisotropic exchange 

As for Cih and D^h symmetry, we first impose the 
invariance of W„e" under tt/2 rotations about the z axis, 
O^ 4 " . For the axial part of H^ 4v , this leads to 



J; 



f,c iv _ 



n,n-\-m 



7i— l,n+m- 



J^-,^4v T 



c,C 4v 



n— l,m+m— 1 



i - Jfs 

= Jc. 



(85) 
(86) 



plus restrictions upon 6%*° +m , 4>n%+ m , and ^tT +m , pre- 
cisely as for C4h and D^h symmetries. 

Now, we impose the mirror plane restrictions of G± v . 
The diagonal mirror plane containing the z axis and the 
line y = -x, — ; — 



iC 4 „ 
~f Z 23 



^ 6 requires O e "" z 12 - - 
By analogy with and £> 4 /j, we choose 6> 12 41 ' = 0, leav- 
ing * 
sinxi2 



A C 4 „ 

ai 2 unrestricted. This leads to x 12 
^ 4 "y, where x^ 4 " = 0f 2 4 " + V^ 4 "- Imposing the xz 



cos xfa 4 ^ 



~c 4 „ 



aC 4 „ 



mirror plane restriction, 0^ 4v x 12 ° — ±x 12 " leads to 
either cosx^ 4 " = or sinx^ 4 " = 0. We then choose 
X12" — 0. A similar procedure is carried out separately 
for the next-nearest-neighbor interactions. 

After checking all of the mirror plane operations of C 4v 
symmetry, we find that the C 4v symmetric anisotropic 
exchange vector group elements may be written as 



:C 4 , 
x n,n+l 

Un,n+1 



~ \ (in - In )X + \ (7n + ln)V, (87) 

2 (7n + 7,7 )i + 2 (Trl - 7n ) V. (88) 
z, (89) 



and the matrix diagonalization parameters satisfy 



} = 



Y C4 M 

'Ai2 ; 



(0,0). 



(90) 



The next-nearest-neighbor vector group elements can 
similarly be shown to satisfy 



V 



Ogfi . The Moriya rules require that 



and 



d 24 4 " are perpendicular to the respective mirror planes 
containing the interacting spin sites. Thus, we write 
d it = d 13 (y ~x)/V2 and = d 24 (x + y)/V2. Then, 
we impose the invariance under tt/2 rotations about the 
z axis, Of 4 ". This operation satisfies Op 41 ' S n = S„ 



and the coordinates transformations O^^x = —y and 
Ci 4 "y = x - This symmetry then forces di 3 = d 24 . 

For the next-nearest-neighbor antisymmetric exchange 
interactions, the relevant mirror planes are the xz and yz 
mirror planes, which bisect the respective near-neighbor 
interaction sites. The Moriya rules dictate that df 2 v and 



»n+i 



d 

the yz plane. We therefore write d 12 



lie in the xz plane, and that d 2 3 4 " and d^ 4 " lie in 



jC 4 „ 
"23 

d C4v 



d 12x x + d 12z z, 
= d 3Ax x + d 34z z, and 



c 

ic, vvu tii^i ciui vvnit; ix 

= d 23y y 

41 — d 4 \ y y + d 4 \ z z. We then impose the cyclic ro- 
tation symmetry 1 4 " , as for the next-nearest-neighbor 
interactions. Invariance of the Hamiltonian under this 
symmetry then requires d\2 X — ~d 23y — —d 34x = d 4 \ y 
and d 12z = d 23z = d 3iz = d 41z . 

Thus, the combined antisymmetric exchange interac- 
tion for tetramers with C 4v symmetry may be written 
as 



n—1 ~ 

■(S n x S n +i) 

2 

+d'J2(y + • ( 5 « x ^+2), (95) 

n=l 

where the three scalar parameters are d 1 — d\ 3 ^/2, 
d z = d\2 Z , and d± = d\ 2x . Thus, for C 4v symmetry, 
there are both site-independent and site-dependent DM 
interactions. The site-independent DM interactions are 
the same as for dh and D±h symmetry. 



:C 4 „ 
X n,n+2 

Vn,n+2 

;C 4 „ 
Z n,n+2 



1 



V2 



{-!$& + 1 n )y), 



V2 



(91) 

(92) 
(93) 



and the matrix diagonalization parameters satisfy 



VII. Si GROUP SYMMETRY 

For S4, the group operations are clockwise and coun- 
terclockwise operation rotations by 7r/4 followed by a re- 
flection in the xy plane, [i^ The structure is pictured in 
Fig. 2(b). These operations are represented by the ma- 
trices 



= (0&",x£") = (O,~). ("4! 



C. Ci v antisymmetric exchange 



(0 s 4 
^1,2 




(96) 



We first consider the next-nearest-neighbor antisym- 



metric exchange interactions, d 



13 



(Si 



d 



c 4 „ 



(S2 x £4). The relevant mirror planes are the diag- 
onal mirror planes containing the z axis and the lines 



A. 54 single-ion anisotropy 



We then construct the vector groups from the site 1 
vectors by successive operations of Of 4 , setting v n+1 = 



10 



;S 4 



v„ , where v n — x n ,y n and z 

-)Si 



~S 4 



The inverse 



operations of 2 are automatically satisfied. With S4 
symmetry, there are no mirror planes to satisfy. How- 
ever, the improper rotation O x l interchanges the x and 
y vector components, complicating the most general vec- 
tor forms. However, we still obtain J^ n = J a 4 n _ 1 = J a 
and J^* = J e „_i — J e , both independent of n, as for 
the three planar symmetries discussed previously. 

Hence, the group vectors for S4 symmetry contain the 
elements 



-Si 
Vn 



|S 4 



where 



COS'i 



/ -cos^f 4 <S+(</>f 4 ) - cosflf 4 sm^f 4 <5,;((/ 
;^f 4 <5-(0f 4 ) - cosflf 4 sin^f 4 <5+(^ 
(_l)n+i S m^f 4 sm6»f 4 



V 



iinV'fM+l^f 4 ) ~ costff 4 cosi/'fM- (. 
sinVf 4 ^(0f 4 ) - costff 4 cos^f 4 <5+ 
(-l)" +1 cosV>f 4 sin6>f 4 



sin vpd n { 
sin0f 4 <5+( 




\ (-l) n+1 cos0f 



K 4 } = 



9f 4 , 



0? 4 >f 4 ) 



(100) 



is unrestricted by the group symmetry, and 8^ {(f)) 
given by Eqs. JjSj! and (J2EJ . 



are 



B. Si symmetric anisotropic exchange 

We now consider the most general symmetric 
anisotropic exchange vector groups for Si symmetry. As 
for Ci V we begin with the near-neighbor axial vector 

Z12 using the general Eq. (|13l) with 9^ iv — > 9f 2 and 
^Civ _^ ^^4_ -y^g cons t ruc t the near-neighbor axial sym- 
metric anisotropic exchange vector group by successive 

operations of Cf 4 , setting £ 2 3 = Ci 4 ^i2i e ^ c - The same 
procedure can be applied to construct the azimuthal sym- 
metric anisotropic exchange vector group. The group 
symmetry imposes no restrictions upon the three near- 



neighbor bond angles 9 



12 1 



fourth operation, Of 4 v^i — 

zSi zSi 

x ,y , and z 



and ijj^, because the 



, ~Si 

for v 



S 4 

v 12 is automatically satisfied 

4 . Thus, the near-neighbor 
symmetry 



anisotropic exchange vector groups with S4 
are fully general, with x n * n+l , y n * n+l , and z n * n+l ob- 

tained from the single-ion x n , y n , and z n vectors in 

-» (A? 4 and 



Eqs. (EE-© by replacing 9** -► 



V 1 ! 4 — * ^12 > respectively. Thus, 



s 4 
12 ■ 



W2 4 } = 



Si 
12 ' ' 



(101) 



is unrestricted by the group symmetry. 

However, the next-nearest-neighbor symmetric 
anisotropic exchange vector groups are different. We 
equal to Eq. JT3J) with flf 4 " -> 0f| and 



first set z^ 

S4 1S4 

6f|, and construct z 2 t from Of 4 z 13 



6f 4 " 



is 



-Si 
'24- 

But we then require Of 4 224 — 231 — ~^13- This forces 
cos^fl = 0, and we take 6*fg = 7r/2. However, 
unrestricted by the group symmetry. 

Performing the same construction of the az- 
imuthal symmetric next-nearest-neighbor exchange vec- 
tor groups, there is a restriction upon ip-ii, which can be 
satisfied by choosing either sin ^3 = or cos -013 = 0. 
We choose sin^fg = 0. We thus obtain 



zSi 

X n.n+2 
Vn.n+2 
Z n,n+2 



-6+(cf>ft)x + 6-(<f>!t)y, 
(~l) n+1 z, 

6-(<t>!$)x + 6t(4>?l)y, 



and 



Si\ 
13 J 



(7r/2,^ 3 4 ,0). 



(102) 
(103) 
(104) 

(105) 



We note that the axial symmetric anisotropic exchange 
vectors are not necessarily collinear with the bond di- 
rections. If all of these anisotropic exchange interac- 
tions arose solely from dipole-dipole interactions, the ax- 
ial anisotropic exchange vectors would lie along the bond 
directions, and the azimuthal anisotropic exchange vec- 
tors would be irrelevant. This would imply cf>f^ = 37r/4 
for the next-nearest-neighbor interactions, and 9f 2 = c / a 
and <f> 12 = 7r for the near-neighbor interactions. How- 
ever, the group symmetry does not require this, and if 
some amount of physical anisotropic exchange (not aris- 
ing from dipolar interactions) were also present, more 
general </>i3,#i2 an d 012 would be allowed, and the az- 
imuthal anisotropic exchange vectors would be relevant, 
with general 1P13 and 1P12 > respectively. 



C. 54 antisymmetric exchange 

We begin with the next-nearest-neighbor antisymmet- 
ric exchange Hamiltonian, which may generally be writ- 
ten as d(g-(Si X S 3 )+df|-(S 2 xS 4 ), where nominally df| 
and dfl are three- vectors. Using the relations Of 4 S n = 

^x = -y,O s 1 \ 
that the components satisfy dfi 



S n+ i, O^x = y, O^y = x, and O^z = -z, we find 



132 



d-24z 



and d°3 X = — dfly There are no mirror planes to satisfy. 
There are therefore two distinct parameters for the next- 
nearest-neighbor antisymmetric exchange. 

For the near-neighbor antisymmetric exchange, we 
nominally set df l 4 n+1 to be general threee- vectors. Then, 
requiring the near-neighbor antisymmetric exchange 
Hamiltonian to be invariant under the improper rotation 
Of 4 , we then find that the component satisfy df 2x = 



0, d% 



d 2 X x i 
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^23y — d 3 l x — dffy, d 12y — d 23x — d 3 ^y — d 4 i x , 
and df| a = -df^ = d%£ z = — df^ 2 . Thus, three parame- 
ters describe the near-neighbor antisymmetric exchange 
for £4 symmetry. 

Combining the near-neighbor and next-nearest- 
neighbor antisymmetric exchange interactions, we find 



= d z z ■ (^(-1)" +1 S„ x S„ 

n=l 

+d ± ■ (Si xS 2 -S 3 x 5 4 ) 
+(d x x z) ■ (s 2 xS 3 -S 4 x Si 

+d' ± ■ (Si x S 3 ) + (d' ± x z) ■ (S 2 x S 4 ), 

(106) 

where d z is a scalar, and d^ and d' ± are two- vectors in 
the xy plane. 

Thus, there are only site-dependent and site- 
independent DM interactions for S4 symmetry. The NN 
DM interactions average to zero, but the NNN DM in- 
teractions generally do not. 



VIII. D 2d GROUP SYMMETRY 

We now consider the interesting case of D 2 d symmetry 
pictured in Fig. 2. In this case, the operations 0® 2d 
A = 1, . . . , 9 are rotations by ir about the x, y, and z axes 
(A = 1,2,3, respectively), and the two diagonal mirror 
planes associated with the principal rotation axis, z.[24| 
The first of these is 



oi 



D 2 




(107) 



and the A = 2, 3 cases have the single +1 in the yy and 
zz positions, respectively. The two mirror planes asso- 
ciated with the principal z axis rotation are the planes 
containing the z axis and the lines y = ±x, respectively, 
Q^l d = Ofg, given by Eq. (HJ. We note, for example, 
that the plane bisecting the sites 1,2 and passing through 
the sites 3,4 is not a mirror plane unless c = a, which we 
assume not to be true for D 3 d symmetry. 

We first try to generate the most general axial vector 

~D 2d 

set satisfying D 2 d symmetry. Taking z l to have the 
general form of Eq. (|13l) . we then require it to be in- 
variant under 0^ 2d , which leaves fi invariant. It is then 



easy to see that we must have 4>® 2d = 37r/4, but Q^ 2d is 

~D 2d 



lD 2 



unrestricted. We then generate the remaining z from 
Of 2d li = I 2 , 2 2d zi = l 4 , and 3 2d zi = z 3 . This 
is consistent with the invariance of the axial part of the 
single-ion Hamiltonian, and also leads to the results that 



We then construct the azimuthal single-ion vectors, be- 
ginning with x 1 ~ d from Eq. (|llfl . using our result that 
0® 2d is arbitrary and 4>f 2d = 3tt/A. Requiring the invari- 
ance of the azimuthal part of the single-ion Hamiltonian 
under the mirror plane represented by Of 2d , we find that 
cosip® 2d — 0, and we choose the solution ip^ 2d = —ir/2, 
so that (j>® 2d + i\)® 2d — 7r /4. Then, we generate the re- 
maining x n 2d from rotations of x 1 ~ d about the x, y, and 
z axes. This of course leads to J® 2d = J e , independent 
of n. 

We then find that the single-ion vector group elements 
for D 2 d symmetry are 



- D 2d I 



-7+cos6»f 2d 



i-D 2d 
Vn 



X " 1 y/2(-\) n am6?* d 



-iy 



V2 



~„D 2d 1 
Z„ 



In \ 








iy 




-7+ sin 


Q D 2d 


7~sin£ 


,D 2d 


»(_l)n+l 


cos 6^ 



(108) 



(109) 



(110) 



V2 



and the associated single-ion parameters are 



K-} = ,-£)■ (in) 



B. Did symmetric anisotropic exchange 

We first consider the near-neighbor axial exchange vec- 
tor between sites 1 and 2, writing £12 in the form of Eq. 
(I13JI . Since we require it to be either odd or even un- 
der rotations about the x axis, 0® 2d Zi 2 — =F^i2, we 
then find that sin</>^ 2d = 0, and we choose (j\ 2 2d = 
with 9® 2 2d arbitrary. We then employ 3 2d zi 2 = 234, 
2 2d Z\ 2 = —234, and Of 2d Zi 2 = —Z41 to generate the 
axial near-neighbor symmetric anisotropic exchange vec- 
tor group. 

We then consider the azimuthal near-neighbor sym- 
metric anisotropic exchange vector x 12 d , using Eq. I jllj l. 
Demanding invariance of the exchange Hamiltonian un- 
der Of 2d , or that 1 2d x l2 — ^x X2 , we then require 
cosV , i2 2d = 0, and choose i\\ 2 2d = ir/2. We then generate 
the remainder of the azimuthal near-neighbor symmet- 
ric anisotropic exchange vector group as above. We thus 
obtain 



zD 2d 
X 



n,n+l 



Ji 



J a -i independent of n. 



( (7^+7rT)cos^ 

(in ~1n) COS °l 

V (-1)" 



(112) 



iD 2l 
'12 



12 



zD 2d 
Hn,n+1 



n,n+l 




\ (-l)" +1 cos6lf 2 2d 
and the associated parameter set is 



K 2 2d } 



,o,- 



(113) 



(114) 



(115) 



d 



jD 2 



1 34y V 



dggz, and d£ 2 ' 



Then, we rotate the near-neighbor antisymmetric ex- 
change Hamiltonian about the x, y, and z axes, and im- 
pose invariance under the diagonal mirror planes. These 
combined symmetry operations result in the following re- 
strictions upon the parameters: d^ d = d^ d — — d^u? — 



a 41x ana "122 — «23z — «: 



jD 2 



-23 Z - "34 W2 = - d 41z d - ThuS > there 

are only two free parameters describing the near-neighbor 
antisymmetric anisotropic exchange interactions. The 
combined antisymmetric exchange Hamiltonian for Z?2d 
symmetry may then be written as 



Next, we construct the next-nearest-neighbor group- 
invariant symmetric anisotropic exchange vector groups. 
Beginning with £13 of the form of Eq. (|13|) . we re- 
quire invariance under the parallel diagonal mirror plane, 
G? 2d z 13 = ±l 13 . This leads to either sin0? 2d = 



or cost 



13 



= 0, and we take 
^d 2 



iD 2l 
7 13 



= 0. Then, 

we similarly require 0^ 2d X\ 3 — ±^13, and require 
sinx^ 2d = cosx^ 2d , choosing x^ 2d = 7r/4, where Xi3 2d = 
4>iz d + V'i3 2d - Thus, the next-nearest-neighbor symmetric 
anisotropic exchange vector group elements are found to 
be 



Xn : n+2 

Un,n+2 = 
Zn,n+2 = 



1 
1 

(_!)n+l 



(■y+x + j-y), 
i. 



and the associated parmaeter set is 

{t4f) = 



qD 2 , 
7 13 



xf 3 2d ) 



(0, ? ). 



(116) 

(117) 
(118) 

(119) 



We note that these vector groups are consistent with 



f;D 2 



n,n+m 



J/, TO and J, 



n,n+m 



Jc,m, as for the previous 



tetramer group symmetries studied. 



C. D2d antisymmetric exchange 



n-/D 2d 
n DM 



= J2( S n X S n+1 ) • ((-l) n+1 d z Z 



-%[(7+ + 7„)* + (7^-7„)y]) 
2 

+d' ± J2l* + (-i)"y] • (s n x s n+2 ),(i20) 

Tl=l 

where d' ± is a scalar and we have set the other two scalars 



d° 2 2d and d A 



d i2y- The DM interactions for D 2d 



symmetry are similar to those for £4 symmetry, except 
that the S4 two- vectors, d± and d±, are just the one- 
vectors, d±y and d'^x. 



IX. T d GROUP SYMMETRY 

Aside from the restriction c = a, T<j is a very high sym- 
metry group. [24] in addition to the identity operation, 
there are 23 other distinct operations. [24^ all of which 
are necessary to include for our purposes. The opera- 
tions 0^ d for A = 1,5 are the same as for D 2 d symmetry, 
rotations by 7r about the x, y, and z axis, and mirror 
planes containing the z axis and the lines y = ±x. There 
are also four additional mirror planes, two of which con- 
tain the y axis and the lines z — ±x (A = 6, 7), and two of 
which contain the x axis and the lines y — ±z (A = 8, 9). 
These may be written as 



We first consider the next-nearest-neighbor antisym- 
metric anisotropic exchange, which couples the spins at 
sites 1,3 and at 2,4, which lie parallel to the xy plane, and 
hence, the roles of the diagonal mirror planes are iden- 
tical to that with C^ v symmetry. Thus, e£^ 2d ||(y — x) 
and d®l d \\(i) + x). Then, we require invariance under 
rotations about the x, y, and z axes. This leads to a sin- 
gle scalar parameter describing the next-nearest-neighbor 
antisymmetric anisotropic exchange interactions. 

The near-neighbor antisymmetric anisotropic exchange 
Hamiltonian is obtained differently. We note that each 
position vector between each of the near-neighbor sites 
has an axis of rotation that bisects it. Thus, d® 2d +1 must 
be orthogonal to the respective axis of rotation. Specif- 
ically, df 2 2d = < 2d y + d?£ d z, A 



23 



d 23x & 



d° 2d z 



n Td 




(121) 



(122) 



In addition, there are four clockwise and four counter- 
clockwise rotations by 27r/3 about the cube diagonals. 
These may be written as 



Cl Td 
^10,12 




(123) 
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(7 Td 

J - y 14,16 




(124) 



and 0^ +1 = [O^j = [O^j for n = 5,6,7,8. 

In addition, there are the six S4 symmetries, which are 
clockwise and counterclockwise rotations about the x, y, 
and z axes. [3 The first two of these are Of| 19 = Of 4 2 . 
The S*4 operations in which the axis of rotation are the 
x and y axes, respectively, are 



(0 Td — 
^20,21 ~ 



^22,23 — 



-1 


\ 







±1 , 


(125) 





Tl J 







±1 \ 







10, 


(126) 


Tl 


0/ 





respectively. 



vector groups, each with six elements. From Eqs. I|114|) 
and i|118[l . the axial vector group {z^ d n+m } is found to 
have the six elements 



v/2 

+(-i)™ +1 



\ ((7n + In )& + iln ~ In )v) 



Om,2 



[x + (-l)" +1 y].(130) 



It is easily verified that this axial vector set is a group 
satisfying all of the Td symmetries. The reason this is 
true is actually pretty obvious: the axial exchange vec- 
tor group consists of equally normalized differences be- 
tween position vectors, and each position vector neces- 
sarily satisfies the group symmetry. However, there are 
no azimuthal anisotropic exchange groups for Td, for the 
same reasons as for the single-ion group vectors. For Td, 
there is therefore only one anisotropic exchange interac- 



tion J, 



f,T d 



ri'n+m 



Jf,i = Jf,2, since J CjTO = for m = 1, 2. 



A. Td single-ion anisotropy 

We first examine the axial single-ion vector for site 
1, making use of Eq. 1|13|) . We require that it remain 
invariant under the mirror planes Ojg s . This forces the 
equations cos 9 — 1 = sin^ 1 d sin0 1 d = — sin 9 1 d cos 4> 1 d , 
which has the solution <fi^ d = 'Sir /4 and 9j d = tan~ 1 (v / 2). 

Thus, Zy = T\. By successive rotations about the x, y, 



and z axes, we generate the Td axial vector set. 
obtain 



z Td 



= - 7 +* + 7,7y + (-i) n+1 i 



We 

(127) 
(128) 



and the associated two-parameter set is 

K d } = (^ d ,^ d )-(tan- 1 (^),3V4). (129) 



Since these axial vectors are identical to the position vec- 
tors, they automatically satisfy all of the remaining sym- 
metry operations. Of course, J Td n = J a is clearly in- 
dependent of n, as it was for the lower symmetry Did 
group. 

We have attempted to construct an azimuthal single- 
ion vector set, but there is no set other than the axial 
single-ion vector set that can satisfy all of the symme- 
try operations. Hence, there are no azimuthal single-ion 
vector groups, and we must infer J Td n = 0. 



C. Td antisymmetric exchanges vanish 

We now consider the possibility that there might be a 
non-vanishing antisymmetric exchange Hamiltonian for 
tetramers with Td symmetry. From the mirror planes, the 
Moriya rules immediately restrict d§ = d\ 2 (z — y)/V2, 
<*23 = d 2 3 (i + x)h/2, dlt = d 3i (z + y)/V2, dg = 
d A1 {z - x)/y/2, d^ = di 3 (x - y)/V2, and dg 4 = 
d2i(x + y)/y/2. Then, rotation by tt about the x, y, 
and z axes (0\ d 2 3) forces ^24 = rfi3, dn ~ d 23 , and 
d\2 = ^34- Then, we require the Hamiltonian to be invar- 
ian under the clockwise rotation by 27r/3 about the cube 
diagonal passing through site 1, 0\ 2 . This operation 



O l2 S2 = S3, 12 Ss 



satisfies 12 x 



x, O^z 
S i7 and O^Si 



y, o^s, 



Si, 

S2 ■ Invari- 

ance of the Hamiltonian then requires <ii2 = c?23 = — di.3. 
This reduces the number of parameters to just one, which 
we take to be d\ 2 . We note that the rotation about the 
cube diagonal O12 leaves this form of TCj^ M invariant, 
without any further restrictions. However, we must also 
have 'H^f M invariant under the six £4 operations, and we 
choose Ofg 

see that setting Oj^Hj) M 
Hence, there are no DM interactions in tetramers exhibit- 
ing T d group symmetry, 



O x 4 as the first example. It is then easy to 
TrJ) M forces di 2 = — di 2 = 0. 



X. THE HAMILTONIAN IN THE MOLECULAR 
REPRESENTATION 



B. Td symmetric anisotropic exchange 



A. The molecular single-ion Hamiltonian 



For Td symmetry, all six of the bond lengths are iden- 
tical, as c/o = 1. Hence, there are at most three local 



To make contact with experiment, we use Eqs. (|39fl - 
gU, §3-©, ED-®, P7| l-iP |l . fTUg|) - ((TT0t . and 
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(|128fl to rewrite Ti 9 f in the molecular (x, y, z) repre- 
sentation, 



H 9 = - 



£(-J#(a«M„ + (-i) n - s 2 n>y ) 

n 



the effective interactions can be large, even if the molec- 
ular structure is nearly Td (that is, if c/a w 1). Once 
this symmetry is broken, the site-independent single-ion 
interactions Jf can become large. 



B. Symmetric anisotropic exchange in the molecular 
representation 



where a, /3 = x,y, z and the /if = {/if} are given by Eqs. 
63, JEU, ifTUOjl . ifTTTl) . and In Eq. (USTJ, 

{A, i?} = A£? + BA and we subtracted an irrelevant con- 
stant. The easiest groups for which to evaluate the pa- 
rameters in T-L 9 si are D 2 d and Td- The site- independent 
interactions in the molecular representation are 



J z 

Jf-(/if 4 ») 



Jf 4h - Ja, (132) 

i (j (3 cos 2 9^ - 1) + 3J e sin 2 

(133) 



^(3cos 2 ^-l) 
-^J e sin 2 Of* cos(2^f 4 ), 



(134) 



- (j a (3 cos 2 flf 2d - 1) + 3 J e sin 2 6»f 2d 



■/J d (Ml 



0. 



(135) 
(136) 



We note that the form of the site-independent interac- 
tions for D 2 d and C 4v symmetry are identical, and can 
be obtained from that for the lowest symmetry, S 4 , by 
setting i/if 4 — > 37r/4. In addition, that for the higher sym- 
metry cases Dih and C 4 h are also obtained from that for 



S 4 by setting 



0. The vanishing site-independent 



single-ion interaction for the highest Td symmetry is ob- 
tained from that of the lowest S 4 symmetry by setting 
both 6>f 4 -> tan-^-y/S) and J e ->■ 0. 

In addition, "H - contains the site-dependent interac- 
tions (-l)"J| !/ (/if) for g = C 4h ,S 4 and K 9 ap (n,^{) for 
g = CihiD^h, C 4v , S 4 , D 2 d and Td, which are listed in the 
Appendix. These site-dependent interactions were gener- 
ated from the site-independent interaction J a for Td, J e 
for C 4 h and DaHi and both J a and J e for C 4v , D 2 d, and 
S 4 in the local coordinates by the respective group oper- 
ations. All of these site-dependent interactions average 
to zero, but they contribute to the eigenstate energies to 
second and higher orders in the interactions J a and J e . 
For C 4v and S 4 , the effective axial site- independent in- 
teractions Jf (/if) arise from a combination of the local 
axial and azimuthal interactions J a and J e , the details of 
the combination depending upon the precise effects of the 
ligand groups near to the ions. We also note that there is 
a substantial difference between Td and the other three 
symmetries. For Td symmetry, the average effective inter- 
action vanishes. However, for C 4v , D 2 d and S 4 symmetry, 



We then use these local exchange coordinates to con- 
struct the group-satisfying anisotropic exchange Hamil- 
tonian in the molecular coordinates. For all four symme- 
tries, there are similar effects. They first lead to renor- 
malizations of the isotropic exchange interactions, mod- 
ifying Hi to 



Hq T = -^S 2 - 7 B-S~ {J ' 9 J9 \ s 2 
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where 



Jg+SJg, 
j'g+SJ' g , 



(137) 

(138) 
(139) 



where J' g is given by Eqs. © and (|10fl . For simplicity of 
presentation, we write 



J, 



m,± 



(140) 



for m — 1,2. Then, the isotropic exchange renormaliza- 
tions may be written as 



6J c ih 

SJ D 2d 

SJ Si 
Sj!p, 



5J D 
1 



0. 



[Jf,i + Jcicos^ia 4 ")], 



-Ji^sin 2 6»f 2 2d , 



sin^6» 



2 nSi 



12 



[J/.1 + J c ,iCos(2^f 2 4 )] 



4 ' 

SJ's 4 
2 



0, 



(141) 
(142) 
(143) 

(144) 

(145) 

(146) 
(147) 

(148) 



They also lead to addition interactions SH 9 e in the molec- 
ular frame. For g — C 4v , S 4 , D 2 d, and Td, these may be 
written as 



2 6-2m 



m—1 n—1 
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+E^ ) * s ,w 3 )(-i)" +1 

x \ Sn,xS n -\-2,y + S n> yS n -\-2,x\ 
n=l a,/3 



We find 



(149) 



where a,(3 = x,y,z, = {0? 2 , (p^ipu}, and fi 9 13 = 
{8f 3 , 4>i 3 , ip 9 3 }. The site-independent anisotropic ex- 
change interactions are found to be 



T C iv r,.C iv 



T^4h T 



(150) 



Ji.T (Mi2 4 ") = 2^/,i+3J c ,icos(2^ 2 4 1: (151) 
= ^(l-3cos 2 ^) 



+ ^Bin a ^co8(2V»S), (152) 
j[J/,i +3Ji,_sin 2 6»f 2 2d ], (153) 



J l,z — A I> 1 ~ z,,J 2,z> 



X 



2,z 



jD^h jD-2d 

J 2,z — J 2,z 



J. 



Si 

2,2 



^2,z" — ^2.+ + J, 



2.+ "t" Jc,2, 



(154) 

(155) 
(156) 



where Jf z (^f 3 ) — Jf z has fully group-restricted ^tf 3 sets 
for all symmetries. The site-dependent anisotropic ex- 
change interactions are listed in the Appendix. 

As we shall see in the following, to first order in the 
single-ion and anisotropic exchange interactions, only the 
site-independent interactions appear. 



d C z lh {n) 


= df 4h (n) = d^"(n) = d z , 


(158) 




= d^(n) = d z (-l) n+1 , 


(159) 




= = L), 


(160) 


d ± 


= d±x, 


(161) 


d ± 


= d±y, 


(162) 


«_L 


= d± = d± :X x + d± :V y, 




«_L 


— «J_ — U > 


(A RA \ 

(164) 




= d' ± y, 


(165) 


d'f™ 


= d' ± x, 


(166) 




= d' ± = d' ±x x + d' ±y y, 


(167) 


d T /{n) 


= d T l = d'f d = 0, 


(168) 



Tetramers with the lowest-symmetry group 54 require 
five parameters to describe the full DM interactions, 
those with either D 2 d or d v symmetry require three 
parameters, those with either C^h or D^h symmetry re- 
quire just one parameter, and tetramers with T c i symme- 
try do not exhibit any DM interactions. We note that 
for C± v , S4, and D 2c t, the NNN DM interactions, while 
site-dependent, do not average to zero. 



XI. EIGENSTATES OF THE FULL 
HAMILTONIAN 

We assume a molecular Hamiltonian of Ti = Hq' t + 
7i a si + STlae- To take proper account of B in 7io, we 
construct our SMM eigenstates in the induction repre- 
sentation by 



C. Antisymmetric exchange Hamiltonian 

In Sees. IV-IX, we already presented the antisymmet- 
ric exchange Hamiltonians in the molecular representa- 
tion. These may be combined into a single equation as 



K B DM = ^(S n xS„ +1 )-(df(n)£ 



/ 




' COS 


dcoscj) 


— sin $ 


sin 8 cos 




•)■ 


COS 


9s'm(j) 


COS(j) 


sin 8 sin 


V 




I - 


sin# 





cos6> 




(169) 

so that B = Bz. A subsequent arbitrary rotation about 
z does not affect the eigenstates. |4j We then set Ti — 1 
and write 



+2(7,7 -7+K 

+ ^(7 + +7„-)(^xdi)) 

+ -J2(Sn * S n+2 ) ■ ~ J+)d'f 

n=l 

+(ri + 7n)(2xdl B )), (157) 

where the relevant parameters are given in Eqs. ( 1591 . 
G3, fTUBTl . (jnOJ, and they all vanish for g = T d . 



s 2 w; m S2i 

s 2 24\r s ^ S2i 
SzW; m S2i 



S24(S24 + 1)|V'I, 1 ^ 24 ), 

mw; m S2i ), 

y/{s — m)(s + m + 1), 



(170) 
(171) 
(172) 
(173) 
(174) 
(175) 



where S& = S$ + iaSy with a = ±. For brevity, we de- 
note v = {s,m, S13, s 24 , {s„}}, and write \v) = IV'f,™' 524 )- 
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From Eqs. JT7U|) to fT75|> . {v'\Hl' r \v) = ££ Q <5„>, where 

E% = -i s ( s + l)- 7 B m 

_ (4- J 9 ) [si3(si3 + 1} + Sm(s24 + 1)tL7g) 

where the J g and are given by Eqs. QJ, i|l(J|) . and 



A. Induction representation 



Hg e to the 



We then transform the molecular TL 9 si 
induction representation, yielding 'H 9 si + 'H 9 ae , and make a 
standard perturbation expansion for the six microscopic 
anisotropy energies {Jj} = {J a , J c ,m, J e ->Jf,m) for m = 
1, 2 small relative to \ J g \, \J'\ and/or 7-B.[4| To do so, the 
contain more interesting physics. Compact expressions 
for these matrix elements are given in the Appendix. 

At arbitrary (0,</)), the first order corrections E 9 ^ = 
{v\H 9 si +'H 9 ae \v) to the eigenstate energies for g = C± v , Si, 
Z?2d, and Td symmetries are 

[3m 2 - S ( s + l)] - - 



a vA — 



J 9 '" \jx 9 ) cos 2 0, (177) 
(178) 



sJ^m = Jl^l)b±-\jlM2)H + K) 



(179) 



where the oi, 6^, and are given in the Appendix, 
and the interactions are given by Eqs. H133|) - (|134[l an d 
(|151fl - <|156|) . We have used the notations 

Jj- 9 = {BLfart^tizM, (180) 

V = i>,Sl3,S24,Sl}, (181) 

which excludes m. We note that the site-dependent in- 
teractions in Eqs. !|131|) an d (|149|l and all of the DM 
interactions in Eq. (|157fl do not contribute to the first- 
order correction to the eigenstate energies. Second or- 
der corrections to the energies include contributions from 
both the site-independent and site-dependent interac- 
tions, and will be presented elsewhere. [3 lj 

For all four g symmetries, E 9 ^ has a form anal- 
ogous to that of the equal-spin dimcr in the ab- 
sence of azimuthal single-ion and anisotropic exchange 
interactions. For these high-symmetry tetramers, to 
first order in the anisotropy interactions, the azimuthal 
single-ion and anisotropic exchange interactions merely 



renormalize the effective respective site-independent ax- 
ial interactions. They do contribute directly to the 
second order eigenstate energies, but in rather compli- 
cated ways, which we will present elsewhere. |3 1| Thus, 
to first order, we only have two effective isotropic ex- 
change and three effective anisotropy interactions, J g , 
J' g , J|(/xf), Jf z (Mi2); an d J'z! which are fixed for a par- 
ticular SMM. Nevertheless, the first-order eigenstate en- 
ergies El = E 9 fi + E 9 'f , given by Eqs. fPToJ) and (|T77|l . 
contain these five effective interaction strengths in ways 
that depend strongly upon the quantum numbers 77 and 
9. These different V, 9 dependencies can be employed 
to provide definitive measures of the five F-independent 
effective isotropic exchange and anisotropy interactions. 
To determine the precise values of the microscopic inter- 
actions, one needs to study the second order eigenstate 
energies in detail. 



B. AFM level crossing inductions 

For AFM tetramers, J g < 0. There will be 2si + 1 
level crossings, provided that the lowest energy state in 
each s manifold does not exhibit level repulsion. Letting 
E f,m ( s i3; s 24, S\) = E 9 + E 9 '^ , these occur when 

E 9 s f 9 («13,*34,*l) = ^, s _l( S 13^24, Sl ). (182) 

There are two types of AFM level-crossings, depending 
upon the sign of J' g — J g , as seen from Eq. (|176Jl . For 
Type I, J g < and J' g — J g > 0, the lowest energy state 
in each s manifold occurs for S13 = S24 = 2s 1; their 
maximum values. Hence, the tetramer effectively acts as 
an AFM dimer with equal spins 2s\. The parameters for 
Type I are given for general s, s\ in the Appendix. For 
the frustrated Type II, J g < and J' g — J g < 0, the level 
crossing occurs for the minimum values of S13 and S24. It 
is easy to see that for even s values, these minima occur 
at S13 = S24 = s/2. For odd s values, the situation is 
doubly degenerate, and the minima occur at 513,524 = 
(s ± l)/2, (s =p l)/2. Hence, this level-crossing behavior 
differs substantially from that of equal spin dimers.|J| In 
the Appendix, we also listed the formulas for the level- 
crossing parameters for Type II. 

For both signs of J g — J' , the sth AFM level-crossing 
induction in the first-order approximation may be writ- 
ten as 



-J g s + (Jg - J'gMJg - J' 9 ) E (c^-) 
- J JM±(a+ + 2b + + a+cos 2 0) 
- J M^[c- + l -{b+ + b-) + c- lC oJ9)] 



J 9 

+ (a 2 +b + a l cos 



(183) 
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where 0(s) is the standard Heaviside step function, E(x) 
is the largest integer in x and 



if = s(2s - l)a£± ,, 



(s-l)(2 S -3)a^± 13iS24 , (184) 
(«-l)a££ u ,.*> ( 185 ) 



a = sa 



^ = "Cm*,^, ( 186 ) 

cr = s (2*-ik-- S24 

_( s _l)( 2s -3KL Ml3 , S24 , (187) 



(188) 



These a^, 6 ± , and cj for j = 1,2 are functions of s, Si 
and the case. For Type II, the functions are different for 
even and odd s. These functions are given for both cases 
with arbitrary s, s± in the Appendix. 

Most notable is that for Type I, the symmetric 
anisotropic exchange interactions combine into a simple 
form, yielding the effective interaction J 9 s (si) given by 



J5f(*i) 



JIM*) , siJL 



4si - 1 



(189) 



Moreover, the symmetric anisotropic exchange contri- 
butions to Eq. I|183|l combine to yield the single term 
— Ci J ( f ff (si)(l — 3cos 2 9)/3 for Type I, where c\ is a func- 
tion of s, si . However, the single ion interactions are more 
complicated, even for Type I, as the ^-independent and 
^-dependent contributions have different s, si dependen- 
cies, except for = 1/2, for which they both vanish. 

One simplification that does occur for Type II AFM 
tetramers is the level crossing that occurs as a re- 
sult of near-neighbor anisotropic exchange interactions 
Jf ^(Mia) • As shown in the Appendix, the contributions 



g,lc(l) _ T g 



^i*(Mi2)/i,*M to jB 9 s ^ s {1) (9) from these 



~fB 

terms are independent of s%, where 



Cl 



^(l + (2s-l)cos 2 
^ Ty (l + (2s -3) cos 2 e), 



s odd 

(190) 
s even. 



Note in particular that for s = 1, f\. z {l,&) — 0, inde- 
pendent of s\. However, the single-ion and next-nearest 
neighbor anisotropic exchange contributions to the level 
crossing inductions depend both upon s, s\ in different 
ways for both odd and even s, as shown in the Appendix. 



si = 1/2 AFM level crossings 

For the simplest case si = 1/2, as in AFM CU4 
tetramers, the single-ion terms in Eq. (|183|l vanish for 
both cases, as for the dimer of equal si = 1/2 spins. |4| 
Using the results given in the Appendix, the expressions 



? 9,lc(l 
1/2, s 

the effective dimer Type I, J' 



for the 7-8^2^'' (0) functions are particularly simple. For 



J g > 0, 



lB{%f{9) = -J 9 + -J c 9 ff (l/2)(l-3cos 2 0),(191) 



1 



G 



jB^\6) = -2J 3 --J c yi/2)(l-3cos 2 0)(192) 



where Jf ff (l/2) is given by Eq. (|189|l with si = 1/2, and 
for the frustrated Type II, J g — J' > 0, 



•^+4^(1 + cos 2 . 



-Jo 



J' 3 



i 9 

J 2,z 



'1 



(193) 



?Vl94) 



Even in this simplest of all tetramer cases, for which 
single-ion interactions do not affect the thermodynam- 
ics, there is still a qualitative different between Type I 
and Type II AFM s\ = 1/2 tetramers. For Type I, there 
is only one effective isotropic interaction, J g and one ef- 
fective anisotropic interaction, J ( f ff (l/2) that affect the 
level crossing. However, for Type II si = 1/2 tetramers, 
the level crossing is different for the near-neighbor and 
next-nearest-neighbor anisotropic exchange interactions, 
and both effective isotropic interactions J g and J' g affect 
the spacing of the level crossings. The only effect of the 
group symmetry is to provide the particular value of J 9 S . 
For example, (1/2) = J/4/8. These expressions also 
show that Type II has a more complex level-crossing in- 
duction variation than does Type I, as the level-crossing 
behavior of Type I is fully described by two parameters 
J g and J 9 S , whereas the level-crossing behavior of Type 
II depends J g , J' g , Jf 2 (/xf 2 ) , and J| 2 separately. On 
the other hand, for this special Si = 1/2 example, the 
9 dependencies of the first and second 7-Bjy2 « are °P~ 
posite in sign for Type I, but have the same sign for 
Type II. In Fig. 4, we illustrated these behaviors for 
J 9 s (l/2)/J g = 0.2 and for the general Type I and for 
Type II also with JlJJg = 0.1 and J g -J' = 0.5| J g \. 



si > 1 AFM level crossings 

For si > 1/2, single-ion anisotropies are allowed, 
and the situation becomes much more complex than for 
si = 1/2. We first consider the simplest si > 1/2 case, 
si = 1, appropriate for AFM M4 tetramers. Exact ex- 
pressions for the s — 1, 2, 3, 4 first-order level-crossing pa- 
rameters are given in the Appendix. In Figs. 5 and 6, we 
plotted the ^-dependence of the first-order level crossing 
induction 7-B? , 1 '(0)/|J fl | for three special AFM cases 
of the general Type I, J' g — J g > 0, and for the particu- 
lar Type II example, J g — J' g = 0.5|J 9 |, respectively. In 
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FIG. 4: Plots of the si = 1/2 first-order level-crossing 
TBf/a^W/l-^l with Jls/Jo = 0-2, J 9 - J g > (solid, Type 
I) and JlJJ g = 0.1, J g - J' g = 0.5|J 9 | (dashed, Type II). 



FIG. 5: Plots of the first-order level-crossing B^ c(1) (0)/|J g | 
for Type I, J g — J' g < and si = 1. Solid curves: Jf (/if )/J g = 
0.2. Dashed curves: J e 9 ff (1)/J 9 = 0.2. 



each curve, we allow only one of the anisotropy interac- 
tions (or effective interactions) to be non-vanishing. In 
Fig. 3, the solid curves are for J|(/ii)/J H = 0.2 for g = 
C\h, £>4/j, C± v , S4, and D 2d . The dashed curves in Fig. 5 
arc for J% s {l)/J g = 0.2 with g = C ih , D 4h , C 4v ,S 4 , D 2d , 
and T d . For g = T d , these dashed curves correspond 
to Jf^/(2AJ g ) = 0.2, a huge axial symmetric exchange 
anisotropy, but still accessible in our first-order pertur- 
bation due to the small coefficient. We note from Fig. 5 
that for Type I, the single-ion and symmetric exchange 
anisotropies lead to opposite ^-dependencies, both hav- 
ing a change in sign just before the second level crossing. 

In Fig. 6, we illustrate the si = 1 level crossings for 
Type II, setting J g — J' g — 0.5 1 J g \. The solid curves are for 

JfUvD/Jg = 0- 2 f° r 9 = Cih,Dih,Ci v ,Si, and D 2 d, as 
in Fig. 5. Except for g — T dl Type II does not exhibit a 
unique effective symmetric anisotropic exchange interac- 
tion, so we plotted the two different symmetric exchange 
anisotropy interactions separately. The dashed and dot- 
ted curves are for Jf Z {nf 2 )/Jg = 0-4 and J 2z /J g = 0.4, 
respectively, which apply for g — C^h, D\h-, Ca v , S4, and 
D 2d . The dash-dotted curves correspond to J/,i/ J g = 0.4 
for g — T d . We note that for each curve, the Type 
II isotropic exchange parameters lead to a larger gap 
between the s — 2 and s = 3 level crossings. The 
sign of the ^-dependencies of the single-ion (solid) curves 
changes between s = 2 and s — 3. The NN symmetric 
anisotropic exchange interactions vanish for s = 1 , but in- 
crease in magnitude with increasing s for s = 2, 3, 4. The 
sign of the ^-dependence of the level crossing due to the 
NNN symmetric anisotropic exchange interactions does 
not change, but its magnitude increases monotonically. 
For T d , there is a sign change in the ^-dependence of be- 
tween the s = 3 and s — 4 level crossings. Thus, Type 
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FIG. 6: Plots of the first-order level-crossing 7-B?^ c(1) (0)/|J g | 
for Type II with J g — J' g = 0.5|J 9 | and si = 1. Solid curves: 
J 9 z {nT)/J g = 0.2. Dashed curves: JL (/i? 2 )M = 0.4, g = 
dh, D4h, C\ v , Si, D2d- Dotted curves: J| z /J g = 0.4 for g = 
Cih, D 4h , C\ v , S4, D 2 d- Dash-dotted curves: J/,i/J 9 = 0.4 for 
9 = T d . 

II AFM si = 1 tctramers have a richer set of first-order 
level-crossing behaviors than do Type I AFM s\ = 1 
tetramers. 

In Figs. 7 and 8, the analogous Type I and Type II 
AFM level crossing inductions are plotted versus 8 for 
si = 3/2 equal-spin tetramers, such as Cr4. The notation 
is the same as in Figs. 5 and 6. For Type I si = 3/2 
AFM tetramers, the single-ion and symmetric anisotropic 
exchange interactions lead to different ^-dependencies of 
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FIG. 7: Plots of the first-order level-crossing ■yB 9 '^\9)/\J g \ 
for Type I, J g — J g < and si = 3/2. Solid curves: 
Jz{l4)/Jg = 0.1- Dashed curves: J e 9 ff (3/2)/ J s = 0.1. 



XII. THE SELF-CONSISTENT HARTREE 
APPROXIMATION 

The self-consistent Hartree approximation, or strong 
exchange limit, provides accurate results for the B 
dependence of the specific heat and magnetization at low 
ksT/\ J\ and jB/\J\ not too small, where k B is Boltz- 
mann's constant. In this approximation, E® = E^ +Ef, 1 
is given by Eqs. I|17(j|) an d l|177|) . respectively. We shall 
present the self-consistent Hartree approximation of four 
measurable quantities in the induction representation. 
We first define the trace in a way sufficiently general as 
to encompass the four group symmetries under consider- 
ation, 



2si 



S13+S24 



E = E E E > 

V S13,S24=0 m = -S S=|si3-S 2 4| 



(195) 



The partition function in the self-consistent Hartree ap- 
proximation may then be written 



Tr^e 



(196) 



CD 




s = 4 



s = 2 



s = 1 



AFM s = 3/2 Type II 



15 



30 



90 



(deg) 



FIG. 8: Plots of the first-order level-crossing jB 9 '^ {6)/\ j g \ 

for Type II with J g - J' g = 0.5|J 9 | and si = 3/2. Solid 

curves: J 3 z {n{)/J g = 0.2. Dashed curves: Jf jS ,(/jf 2 ) / J g = 0.4, 

g = Cah, Dak, Cav, Sa, D2d- Dotted curves: J| z /J 9 = 0.4 for 

g = Cah, Dah, Cav, Sa, D2d- Dash-dotted curves: Jf,i/J g = 
0.4 for g = T d . 



the level-crossing inductions, each with a change in sign 
in the 9 dependence at about the second level crossings, 
as seen in Fig. 7. For Type II s± = 3/2 AFM tetramers, 
the sign changes appear between the second and third 
level crossings, as shown in Fig. 8. 



where (3 = l/(k B T), and E® = E 9 ufi + E 9 v f is given by 
Eqs. (|176|l and l|177(l . In this compact notation, the self- 
consistent Hartree magnetization Mg X \B,T) and spe- 
.(i) , 



cific heat C^^B, T) are given by 



mW(B,T) = jTr^me-^/Z, 



(i) 
g > 



(197) 



C£1(B,T) 
k B (3 2 



Tx,(^e-^)n 



(198) 



We note that there are strong differences between the 
low-T behavior of FM and AFM tetramers. We assume 



\Jg\> 



\J> 



J g \. For FM tetramers with J g > 0, the low- 



T thermodynamic behavior is dominated by the s = 4si, 
m = — 4si state, leading to 

MW(B,T) t Z jBB isi (0-yB), (199) 

where Bs(x) is the Brillouin function. The universal- 
ity of this function renders thermodynamic studies use- 
less for the determination of the microscopic parameters. 
For AFM tetramers with J g < 0, however, there will 
be interesting level-crossing effects, which can be em- 
ployed to measure the microscopic interaction parame- 
ters, as discussed in detail in Sec. V. As for dimers, 
C V (B,T) for AFM tetramers at sufficiently low T ex- 
hibit 2s i central minima at the level-crossing inductions 
5f; lc s (6>) « B 9 s ^s {1) {9) that vanish as T -► 0, equally sur- 
rounded by peaks of equal height. 4] As for the magne- 
tization, Cv(B,T) for FM tetramers at low T reduces 
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to that of a monomer with spin 4si, yielding a rather 
uninteresting Schottky anomaly. 

However, the microscopic nature of FM tetramers can 
be better probed either by EPR or INS techniques. The 
self-consistent Hartree EPR absorption Qx-a a^, w ) f° r 
clockwise (er = 1) or counterclockwise (a = —1) circu- 
larly polarized oscillatory fields normal to B is 



where M v v * 



7 



YjTr^Tr^e 



-0El 



\M V 



x [&[E° - El, + u) - 5{El, - El + w )] , 

(200) 



Urn' ,m+cr 



y^„/. The strong resonant inductions for appear at 



and Tiv = 



7 B*f) = ± w+ ^f^(l-3cos^)jr(^), 

(201) 

where J%' v {ji 9 ) is given by Eq. I|178|) . We note that 
J| ,!/ (/i 9 ) contains the three effective microscopic inter- 
actions, J|(/if), J% 2(^12) ' anc ^ ^2,2, multiplied by the 
constants a— , — c^, and — o^/2. Thus, for an experimen- 
tal determination of all three microscopic parameters, it 
is necessary to perform EPR experiments in three differ- 
ent quantum state manifolds. At low T, the Boltzmann 
factor favors the ground state manifold, so to obtain the 
desired measurements of three manifolds, one needs to 
access the excited s — 4s 1 — 1 states, as well. 

As shown in the Appendix, for a FM tetramer with 
J g > in the lowest energy manifold of states, nominally 
(s, S13, S24) = (4si, 2si, 2si), has the highest population 
at very low T. In this ground state manifold with — 4si < 
m < 4s 1, the constants reduce to 



s- 1 T 1 
2(2s- 1)' 

s - 1 
2(2s- 1)' 



(202) 
(203) 



where s = 4si. We note, of course, that when the single- 
ion and anisotropic exchange interactions are sufficiently 
large, s is not really a good quantum number, as the 
actual ground state, while mostly the nominal s = 4s 1; 
contains a mixture of s — 4si — 2 states, to leading order 
in those interactions. Q Additional mixtures are likely to 
be relevant for the excited states, as well. 

However, assuming for simplicity that we can safely 
neglect such mixing, and take (s,Si3,S24,) to be a good 
quantum number set, there is still some ambiguity in 
the determination of the order of the lowest manifolds of 
excited states. From Eq. (|176fl . it is evident that if J g > 
and J' g = J g > 0, the states (s,si3,S24 = (s,2si,2si) 
are energetically favored for each s value. We denote 
this a Type I FM tetramer. However, for a Type II FM 



tetramer with J g > and Jg' — J g < 0, the states with 
minimal allowed values of s 13,324 are favored. Thus, for 
even s, these are (s, s/2, s/2), and for odd s are (s, (s + 
l)/2,(*-l)/2). 

Thus, there are two inequivalent lowest energy ex- 
cited state manifolds with s — 4si — 1 and — 4s 1 < 
m < 4si — 1 of excited states with s = 4si — 1. For 



\Jg-J' n \ < J, 



„ g , these are ordinarily both accessible. One 



of these is non-degenerate (at fixed m) with (s, S13, S24) = 
(4s 1 — 1, 2si, 2si) and the other is doubly degenerate with 
(s, S13, s 24 ) = (4si - 1, 2si - 1, 2si), (4si - 1, 2si, 2si - 1). 
In Type I FM tetramers, the (4si — l,2si,2si) states 
have lower energy than the next excited manifold with 
(4si - l,2si - l,2si),(4si - l,2si,2si - 1) states. In 
Type II FM tetramers, the reverse order occurs. From 
the results presented in the Appendix, for the (4s 1 — 
1, 2si — 1, 2si), (4si — 1, 2si, 2si — 1) manifold of states 
with — 4s 1 + 1 < m < 4si — 1, the constants are 



1) T (2s - 1) 



2s(2s - 1) 
(s + l)(s-l) 2 



2s 2 (2s - 1) 



(204) 
(205) 



where s — 4si — 1. For the (4si — 1, 2si, 2si) manifold of 
states, also with — 4si + 1 < m < 4si — 1, the constants 
are 



(s-l)Q»Tl) 
s(2s-l) : 
(s-2) 



2(2s- 1) 



(206) 
(207) 



where again s = 4si — 1. We note that these three sets 
of constants are all distinct. Thus, either by operat- 
ing at higher temperature or by some special pumping 
technique, it is possible to probe these low-lying state 
manifolds, and thereby obtain a full measurement of the 
first-order microscopic parameters. 

The Hartree INS cross-section S^^B, q, u) is 

^-^'Hv\Sl A \v'){v'\S n3 \v) 



xe 



xSicu + El-El,), 



(208) 



, COS( 



>b,q, qy 



where a,/3 = x,y,z, qz = sin( 
sin6*b 9 sin0b ig , and gj = cos^b,g, @b,q and <fib,q describe 
the relative orientations of B and q,Q the r n and 
{v'\S nt a\v) are given by Eqs. QJ, l|255|l . and (|256|l re- 
spectively, and the scalar q ■ (r„ — r n / ) is invariant under 
the rotation, Eq. 1169( 1 . After some algebra, we rewrite 
Sf'(q,u) as 



El -El,) 
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x (sin 2 6 h>q L v y(q) 



(209) 



where the Hartree functions L v y(q) and M v . v i(q) are 
given in the Appendix. They are independent of B. Since 
E 9 is well-behaved as B — > 0, Eq. 1)209(1 is accurate for 
all B. 

As for the dimer,0 additional EPR and INS transi- 
tions with amplitudes higher order in the anisotropy pa- 
rameters Jf (/if ), Jf Z {H\ 2 ), and J 2 z relative to J g are ob- 
tained in the extended Hartree approximation, but will 
be presented elsewhere for brevity. 



pro? 

m 



XIII. DISCUSSION 

We note that the phenomenological total spin 
anisotropy model widely used in fitting experimental 
data on SMM's is 



Tip — Jb^z Jd(Sx Sy 



(210) 



containing axial and azimuthal contributions, 
respectively, 111 sometimes with additional quartic 
terms, [lil l36| The terms are generally defined relative 
to the total spin principal axes, which for equal spin, 
high symmetry systems are the molecular axis vectors. 
It is easy to evaluate Ef, = (v\H p \v) in the induction 
representation. One obtains Eq. ((177(1 with 



jf 
Jd 



Jb, 

0, 

0. 



(211) 
(212) 
(213) 



A non- vanishing Jd would lead to a term in Eq. 11771) 
proportional to sin 2 9 cos(20), as for the dimer,|4j which 
does not arise in the first-order calculation for the high- 
spin tetramers under consideration, based upon the mi- 
croscopic parameters alone. Hence, Jb alone describes 
the 9 and m dependencies of Ef, 1 correctly. 

In addition, for the lowest energy states of FM 
tetramers, we have s = 4si, S13 = S24 = 2si, so 17 
is restricted to a single set of values. In this high-spin 
case, Jb can provide a correct phenomenology of the 
s = 4si ground state energy manifold. However, if ap- 
plied to the two low-lying excited state manifolds with 
s — 4si — 1, for instance, one would infer two different Jb 
values from that obtained in the ground state manifold. 
For AFM tetramers, the phenomenological model also 
correctly provides a vanishing first-order correction to the 
s = manifold of states with s i3 = s 2 4 = 0, 1, . . . , 2s\. 
However, it also has problems describing the first excited 
manifold of AFM states with s = 1, because it leads to 
the unphysical choice of SJ^' w (fb) = (or a constant, 
independent of the quantum numbers V). Hence, the 
phenomenological model works best in describing only a 



single manifold of states with fixed (s, S13, S24). This is 
more restrictive than the usual assumption of its appli- 
cability to a fixed set of s states. [la kSa 1371 138| 

We note that for the FM Fe 4 SMM with D 3 
symmetry, Fe4(thme)2(dpm)g, where Hsthmc 
is 1,1, 1-tris (hydroxymethyl)ethane and Hdpm is 
dipivalovlmethane. |36l l37| the D3 symmetry also pre- 
cludes the Jd term. Nevertheless, in fits to INS data, 
it was necessary to assume Jd =/= 0, 36, 38] to obtain 
the appropriate anticrossing gaps, so that either the 
powdered sample did not have pure D3 symmetry, 
or the phenomenological model they used, Eq. I(210() 
plus two quartic terms obeying D3 symmetry, was not 
appropriate. Since single-ion interactions appeared to 
be important, [lil the total spin might not have been 
a well-defined quantum number, as in at least one Fe2 
dimer and in Fes-0, IE 0> HI In addition, second order 
effects might also provide a possible explanation.^ 
To investigate this symmetry, it would be necessary 
to evaluate the single-ion matrix elements for that 
symmetry in compact form.j^lj 

For simplicity, we have neglected biquadratic exchange 
and higher order anisotropic interactions. Of these, those 
of the simplest form are the isotropic biquadratic ex- 
change interactions, which for the symmetries under con- 
sideration take the form 



V 9 

n bq 



where J9 



4 

-Ji, g ySS-n 

n=l 



2 

71=1 



'n+2 



) 2 . 



(214) 



V 

J 2.T, 



2,9 7^ Ji,g for 9 = C^h,Dih,C iv ,Si, and D 2 d, 
= J2,T d - The J^ g term is diagonal in the the 
\v) representation, but the remaining term proportional 
to Ji. g is not. However, it is important to note that 
these isotropic interactions are rotationally invariant, so 
they are independent of 9 in the induction representa- 
tion. Hence, they modify the positions but not the 9- 
dependencies of the AFM level-crossing inductions. Al- 
though the diagonal J 2 term does not affect the EPR 
transitions, the non-diagonal J 2 , g term will modify the 
positions but not the ^-dependencies of the EPR reso- 
nant inductions. 

We note that our formulation of the single-ion ma- 
trix elements in terms of a pair of dimers is appli- 
cable to low-symmetry systems such as M01 20.sn (^2- 
OH)ioH2{Ni(H 2 0) 3 } 4 , abbreviated as {NLMo^}.]^ 
systems such as Ni 4 tetramers obtained from salts of 
[Ni 4 (H20)2(PW 9 34)2] 10 " with C 2v symmetry^ and 
the unequal-spin systems Mn^Mn^ 77 and Nin 7 Mno J7 . |4ll 
li^ In the first system with C\ v symmetry, one would 
expect many more single-ion, symmetric anisotropic ex- 
change, and DM interactions, making definitive fits to 
the existing powder magnetization data problematic. [39^ 
However, to improve the fits to the Cs v or D3 symmetry 
systems, such as the Fe 4 compound Fe 4 (thme)2(dpm)6 
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and the Cr^Nig 7 tetramer with an s = 9/2 ground 
state, |3^, 0| would require a reformulation of the single- 
ion matrix elements as a trimer plus a monomer. |3l| 



SUMMARY 

We presented a microscopic theory of high-symmetry 
single molecule magnets, including a compact form for 
the exact single-spin matrix elements for four general 
spins. We used the local axial and azimuthal vector 
groups to construct the invariant single-ion and symmet- 
ric anisotropic exchange Hamiltonians, and the molecular 
representation to obtain the Dzyaloshinskii-Moriya inter- 
actions, for equal-spin si tetramers with site point group 
symmetries T^, D^h , D2d, &4, C^h, or C^ v . Each vector 
group introduces site-dependent molecular single-ion and 
anisotropic exchange interactions. Assuming weak effec- 
tive site-independent single-ion and exchange anisotropy 
interactions, we evaluated the first-order corrections to 
the eigenstate energies, and provided analytic results and 
illustrations of the antiferromagnetic level-crossing in- 
ductions. We also provided Hartree expressions for the 
magnetization, specific heat, EPR absorption, and INS 
cross-section, which are accurate at low temperatures and 
arbitrary magnetic fields. For ferromagnetic tetramers, 
we provided a procedure for a precise EPR determination 
of the microscopic anisotropy parameters. Our procedure 
is extendable to more general systems. 

We thank N. S. Dalai for helpful comments. This work 
was supported by the NSF under contract NER-0304665. 



APPENDIX 
1. Site-dependent single-ion interactions 

For g = D^hi the only non- vanishing site-dependent 
single-ion interaction is 



i^(ri, M ?«) = (-l)" +1 J e 



(215) 



For g = C^h, the two non- vanishing site-dependent 
single-ion interactions are 

J* v 4h (p? 4h ) = J e cos(2 X fn (216) 
K°f{n,iif") = (-l)" +1 J e sin(2 x f-). (217) 

For g — C/t v , the non- vanishing site-dependent single- ion 
interactions in Eq. (|131[) are 



/1 \n+l 



J a sin 2 6^ 
Je(l-r-COS 2 0f 4 *)V (218) 



Til 



2V2 



sin(20f^)(J a - J e ),(219) 



uin(20f*«)(J o + J e ).(220) 



For S4, the single-ion site-dependent interactions are 

4 S >f 4 ) = Ji(Mf 4 )cos(2^) 

+J 2 (/if 4 )sin(20f 4 ), (221) 
K^(n,^) = (-l)"(j 1 ( A1 f 4 )sin(2^) 

-J 2 (/zf)cos(20f 4 )), (222) 

K%(n,tf*) = (-l)" +1 (j 3 (/if 4 )5-(^) 

-J 4 (^)<5+(^f 4 )), (223) 

+J 4 (/if 4 )^(0f 4 )), (224) 



Si \ 



J a sin 2 flf 4 



-J e (l + cos 2 6»f 4 )cos(2V'f 4 )), (225) 
J 2 {l4 4 ) = J e cos6»f 4 sin(2V'f 4 ), (226) 



Mvf 4 ) = ^sin(20f 4 ) [J a + J e cos(2^f 4 )] ,(227) 



J 4 (^f 4 ) = J e sin6»f 4 sin(2^f 4 ) 



3a , 



(228) 



For Z?2d, the non- vanishing site-dependent single-ion 
interactions are 



■(j Q sin 2 C 
+J e (l + cos 2 0f 2d )), (229) 

: ' n -.(J a -J e )sm(29? 2d ), (230) 



2V2 



2^2 



(J„- J e )sin(20? M ). (231) 



The site-dependent interactions for symmetry are eas- 
ily obtained from those Eqs. i|229[l - i|231fl by setting 
6>f 2d -> tan- 4 (V2) and J e -► 0. 



2. Site-dependent symmetric anisotropic exchange 
interactions 

For g = Cih, the non- vanishing site-dependent sym- 
metric anisotropic exchange interactions in Eq. I|149|) 
are 

Jm 4 ,xy{P'l 4 m+l) = Jc,m cos (2Xl,m+l) (232) 

for m — 1,2. For D^, the non- vanishing site-dependent 
symmetric anisotropic exchange interactions are 



(233) 
(234) 



For C^v, the non- vanishing site-dependent anisotropic 
exchange interactions are 



J l,xy\l Jj 12 ) 



\[Jf,i- Jc,icos(2^)], (235) 
J 2 ,_, (236) 



M~ %)Jc,i sin(2^)(237) 



M+%)J c ,ism(2^). (238) 



Again, the more interesting group is g = Si. We find 



-J 1 ( A1 f 2 4 )cos(20f 2 4 ) 



Jty(rfl) = -J 2 ,-sm(2^ 3 4 ), 



J Cl _ cos(20§), 
(-l)" +1 fJi(/if 2 4 )sin(20 



+J 2 (/4 4 )cos(20f 2 i ) 



(-l)"(j 3 (Mf 2 4 )^(0f 2 4 ) 
-^4(a4 4 )<5+(0?2 4 ) 
= (-l)"(^(/xf 2 4 )C(^ 



+.W 2 4 )<5. 



(239) 
(240) 
(241) 

(242) 
(243) 
(244) 



where 



^i(/4 4 ) - 



^2 qSa 



- J c ,i(l + cos 2 6>f 2 4 ) cos(2V>f 2 4 )) (245) 
J C)1 cos 6>f 4 sin(2^f|), (246) 
\ sin(2^)[J / , 1 + J c ,x cos(2V§)], (247) 



J 4 (m 54 



12 , 



J c .i sin 6»f 2 4 sin(2^f 2 4 ), (248) 

and the <5^(</>) are given by Eqs. (|3"4*j) and (|3~5)) . and listed 
in Table II. 

For g = -D 2( 2, the non- vanishing site-dependent 
anisotropic exchange interactions are 



Ji,+ + Ji- cos 2 6>f 2 2d , 

— Jc.2, 



(249) 
(250) 



-(7++7^)^Psin(2<-), 



(7 + -7 - ) — 

V in in / a 



sin(20? 2 M ). 



(251) 



(252) 



For g = Td, the non- vanishing site-dependent anisotropic 
exchange interactions reduce to 

xy 4 ' 



J l,x 



(253) 



K T d _ _JfA 

^2,xy - g ' 



23 
(254) 



3. Compact single-ion matrix elements 

By using the Schwinger boson technique ol represent- 
ing a spin by two non- interacting bosons, and check- 
ing our results using the standard Clebsch-Gordan al- 
gebra with the assistance ol symbolic manipulation soft- 
ware, we find the single-spin matrix elements with gen- 
eral {s n } = (si, s 2 , s 3 , s 4 ) to be 



(v'\S n ,- z \v) = 



^ s 
-^{s n }.s,n 

Sl3,s' 13 ,S24,s' 24 



^{s n },s,n 

S 13' S 13 ,S 24' S 24 



Si ,S3,fi 



^Si 3 ,S 2 4 
-51,-53,^ 



{s„, },s,n 

S 13- S 13' S 24' S 24 

{s n }. — s — l.n 

S 13' S 13 ,S 24' S 24 



a 13' a 13' ^24^24 



(255) 



(v'\SnM") = s m ,, m+ *[A* m 5 s ,, s r { s s " } s f< n s s , 

» S 13' S 13' a 24' a 24 



c 7-j cr , rn a {s re },— s— l.n 

-0 S > s+l JJ _ s _iL± / / 

^13' 13' 24^24 



^ |/;: - A !- M,» ) C256) 

a 13' s 13' 24'*24 7 



\/s 2 — TO 2 , 



s 24 ,s,n, , , _ 



(257) 

= cry (s — <rm)(s — am — 1), (258) 

— <5s 2 4^24 £ ™ C t s 2 1 4 s 3 ( s 13i s 13) 

+^ 3>Sl3 e+<^"( S 24,44), (259) 
= ^s 24 ,s 24 e n /3| 4 4 s3'™(s 13 , s' 13 ) 



+^ 3 , Sl3 e+fe s '> 24 ,4 4 ), (260) 

^(1 + fs,s 13 ,s 24 )£ s ' 13 ,s 13 



'»» ^Sl.SS.S °^13' S 13- 

I p«13+l-S24 X , 

Si,S 3 ,S U S 13 ,S13 + 1 



(261) 



-«7. 



/n \^Si,S3,S%3,S 13 - 

I y^fSi 3 + l,S 24 r 
+^ S1>S3 ,- S "< 3 , Sl3 + l 

/1S24 /IS-Sia 
'/S13,ai,s 3 ^ 1 - S + Sl 3^ 1 - S24 

4s(s + l) 

n A S24 A S24 

'/S13 |S1 ,83 ^8+^13 ^-8+813-1 

4sV4s 2 - 1 



A x+z /i y 



y/z^Az 2 - 1)' 



(262) 

(263) 

(264) 
(265) 
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x(x + 1) - y (y + 1 

= / i r\ » 266 ) 

z(z + l) 

where 7+ and A™ are given by Eqs. J5J and l|175(l . The 

prefactors m, A* m , Cf, C™ s _ u Df' m , and lf^_ x are 
consequences of the Wigner-Eckart theorem for a vec- 
tor operator. [24| The challenge was to obtain the coef- 
ficients r^ s "b s '™ „/ and A^ s "b s '™ , . Their hierar- 

S 13' S 13' 24' S 24 s 13' 13' 24' S 24 

chical structure based upon the unequal-spin dimer sug- 
gests that analogous coefficients with n > 4 may be 
obtainable. |4| Details will be presented elsewhere. 



4. First-order eigenstate energy constants 

The constants appearing in the first-order eigenstate 
energies (|177|l are 



Ui±e 

4±2(E[(^ +1 

V=±l 

E/ ( -,Sl3 + (l+(T)/2,S24 



2 _ 2 \ 

"s,Sl3,S 2 4 "»+l|«13>S24 J 



)/2,«24 



)/2 



£ = 



<t'=±1 



,0 \ ( r si3+(i+<y)/2, S24 



<r' = ±l 

+ (si3^S2 4 )^), (267) 

i £ (2 S + l + ( T') 2 f^ + (l + CT ')/2, Sl 3,,24 
I ^Sl3 + (,H-°'j/ z . s 24 ' 2 

\^ Lj Si,s 1 ,(T(t's+o-(1+(t')/2 / 

+ (si 3 ^ s 24 )l ), (268) 



where the F^ 2 *, G a a ™£™, Vz,x, y , and £3^, are given 
by Eqs. respectively. 



5. Type I First-order AFM level-crossing constants 

The Type I relevant both for J' g - J g > AFM level- 
crossing inductions some low energy FM manifold states 
(with J' g > J g > 0) is S13 = S24 = 2s 1 with s arbitrary. 
Since all spins have si, we let a^ Sl3iS24 = a^, etc. For 



Type I, we have 



s,2si ,2si 



(269) 



1.8!, ± 

u s,2si,2si 



s(s+l) 



2(2s + 3)(2s- 1) 
+ [2s(s + 1) - l][8si(2si + 1) - s(s + 1)] 

± 4s7^I ( 2[16s * + s{s + 1)][s(s + 1} " 1] 

-8si[2s(s + f) - 1]) , (270) 



^s,2si ,2si 



3[a(a + !)-!]- 8ai (2ai + l) 
2(2s-l)(2s + 3) 



(271) 



From these expressions, we may evaluate the Type I 
first-order level-crossing inductions for AFM tetramers. 
From the definitions of the level-crossing constants in 
Eqs. 1)184(1 - H188|) . we rewrite them to explicitly indicate 
the Si, s and type dependencies, and for Type I, it is easy 
to show that 

*?3*W = ^"W( 1= F4^)' (272) 
2s[8s x (2si + 1) + 4s 4 - 10s 2 + 3] 
(4s 2 - l)(4s 2 -9) 

X ( 1:F 4^l)' (273) 
3[4s 3 + 5s 2 - 3s - 3 - 8si(2si + 1)] 



2(2s + l)(2s + 3) 



(274) 



C2o(s) + (4s 2 - 4s + 3)8si(2si + 1) 
2(4s 2 -l)(4s 2 -9) 



c 20 (s) = 3(4s 4 -9s 2 -s + 3), 



(275) 
(276) 



for j — 1,2. For si = 1/2, it is easy to see that 
a l,j 



(s) = &y 2,+ (s) = for s = 1, 2, as expected. 



Moreover, it is easy to show that for Type I, 



(277) 



a~ 2+ b- = -|ar = -^1)^.(278) 

This implies that for Type I, the axial near-neighbor 
and next-nearest-neighbor axial anisotropic exchange in- 
teractions may be combined to yield an effective axial 
anisotropic exchange interaction given by Eq. (|189|l . We 
emphasize that this is an exact expression for the effective 
interactions affecting the first-order level-crossing induc- 
tions, and is independent of the level-crossing number s, 
provided that Type I applies. 

For the single-ion contributions to the level-crossing 
inductions, no such simple relation can be found. We 
note that 



a+ + b+ = ~af, 



(279) 



but the overall quantity a t + 26 + + cos 2 9 in Eq. I|183|) 
contains the extra quantity b + , which depends upon s, s\. 



6. Type II First-order AFM level-crossing constants 

For AFM tetramers with J g < 0, J' g - J g < 0, Type II, 
there are two classes of minimum energy configurations 
for each s manifold, depending upon whether s is even or 
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odd. For even s, the relevant states have S13 = S24 = s/2, 
and for s odd, they are S13, S24 = ( s ±l)/2, (s=Fl)/2. This 
type is also relevant for FM tetramers with J g > J' g > 0, 
especially for the first excited manifold of states with 
s = 4si — 1. For even s the relevant parameters are 



(s l)(2s 3)o g ]j 1 ^ g _ 1 ^ 2 ; ( s _i)/ 2 



U s,s/2,s/2 

his, Si) 

fj>i,± 
s,s/2,s/2 

/2(s,Sl) 



1 



IT 



2(2s-l)L 2(s + 3) 

16si(si + 1) - s 2 - 2s + 6, 

' r s 2 ± /2 0, s l) 



2(2s - 1) 



s + 3 

16si(si + l)(s 2 + 2s-l) 
-s(s 3 + 4s 2 + s-4), 
s- 1 



*.*/2,«/2 2(2s-l)' 
For odd s, the relevant parameters are 



(280) 
(281) 
(282) 

(283) 
(284) 



si,± 



1 



*«,( 8 +i)/2,(«-i)/2 ~ 2s(2a-l) 



s 2 - s + 1 



/3(s,Si) 



/3(s,Sl) 



D s,(s+l)/2,(s-l)/2 



T 2(s + 2)(s + 4). 
16si(s! + l)(.s 2 + 3s-l) 
-s 4 -5s 3 + lls-ll, 



2(2s- 1) 

± 



s 2 -l 



(285) 
(286) 

(287) 



/4(s,Sl) 



s,(s+l)/2,(s-l)/2 



2s(2s- 1)' 

(* + !)(* -I) 2 
c «,(s+i)/2,( s -i)/2 2s 2 (2s-l) ■ 



"(s + 2)(s + 4) 
16si( Sl + l)(s 3 + 5s 2 +4s-3) 
-(s + l)(s 4 + 6s 3 + 7s 2 - 3s + 1), 

„2 „ , 1 

(288) 
(289) 



From these expressions, we obtain the level-crossing in- 
ductions for the AFM type J g < and J' g - J g < 0. 
For s even, we have 



s(2s - i)<;^ 2iS/2 

-(s-l)(2 S -3K^ s/2;(s _ 2)/2 , (290) 



«i,± 

l ||o,2 



sa 



Sl,± 

s,s/2,s/2 







\ Sl,± 

" L ) a s-l,s/2,(s-2)/2' 


(291) 


bTi^s) 


= 


- h su± 
s/2 u s-l,s/2,(s-2)/2' 


(292) 


c nti( 8 ) 


= s(2s- 


" 1 ) C fs/2,s/2 





-( s -l)( 2s -3)c±_ 1;S/2;(s 



-2)/2 



-He, 2 



(s) 



"s, s/2, s/2 



(s l) c f-l, s /2,(s-2)/2- 

For s odd, we have 



, (293) 



(294) 



81,± 

l IIo 



\(s) — s(2s l)a s ^(^!|_i)/2 7 ( s -i)/2 



«i,± 

l IIo,2 



usi,± 
J IIo 



(s) 



(s) 



Si,± 

Sa s,(s+l)/2,(s-l)/2 

~{S - l) a s ^l ; ( s _l)/ 2 ,( s -l)/2' 



(295) 
(296) 



D s,( S +l)/2,( S -l)/2 



-6 



Sl,± 



-l,( s -l)/2,( s -l)/2> 

(297) 



Si,± 

c J/o 



1 (s) — s(2s l) C Sj ( s+ l)/ 2 ,( s -l)/2 

— (s - l)(2s - 3)c^_ 1 ; ( s _i)/ 2l ( s _i)/2' 



= sc 



s,(s+l)/2,( s -l)/2 



(s l) C s _l,( s _l)/ 2l ( s _l)/2- 



(298) 
(299) 



From these expressions, we may obtain the Type II 
AFM level-crossing induction parameters. For even s, 
we find 



a c 10 {s) = 



2s- 3 ^ 48s 1 (s 1 + l) + a c 10 (s) 



4(s + l)(s + 3) 



-2s 3 - 5s 2 + 6s + 18, 



(301) 



(s) = 



1 



2s 2 - 6s + 3 



2(2s-l)(2s-3) 
16si(si + l)(2s 2 - 4s + 3) + a| (s) 



± 



2(« + !)(« + 3) 



a2o( s ) 

J IIe 



2s 4 + 2s 3 - 9s 2 - 18s + 18, 
s 

1 



(2s- l)(2s-3) 
16«i(si + l)(s- 2) + &§(«) 



± 



b C o(s) 



4s 4 - 7s 3 + 20s 2 + 14s - 18, 



-He 



2(s + l)(s + 3) 
7, 

s(2s- 3) 

s(2s 2 - 4s + 1) 



^ II ^> 2(s-l)(2s-l)(2s-3)' 
Combining a 7 }o^(s) and b^ 1 /^ (s) , we find 



l IIe,2 + °//c 



d° (s) 



1 16si(ai + l)-dg(a) 

2 4(s + l)(s + 3) : 
2s 3 + 7s 2 + 2s - 6. 



(302) 
(303) 



(304) 
(305) 
(306) 

(307) 



(308) 
(309) 



We note that this expression differs substantially from 
that for except for Si = 1/2 and s = 2. Simi- 

larly, it is elementary to combine cj le 2 (s) and [i> 77 e + (s) + 
6&-W1/2. We find 



W( S ) + K 6 ^ +(S) + ^" (S) ) 



2(5-1) 



(310) 



This simple expression differs from that for c u 1 (s) by 
the factor 2s — 3. However, at s = 2, the only even s 
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value for si = 1/2, they are equivalent. In addition, as 
for Type I, there is no simple relation between the single- 
ion parameters a^j^l-s) + 2&^' c + (s) and <Xjjei( s ) - 

For odd s, the Type II level-crossing induction param- 
eters are 



si , ± / \ 



a° w (s) 



2s 



2 

-2s 3 



1 _^ 48si(si + 1) + a° 10 (s) 



4(s + 2)(s + 4) 



lis 2 



± 



10s + 17, 

2s 2 - 2s - 1 
16si(si + l)(2s 2 + l) + a§ (s) 



2(2s- l)(2s-3) 



a%(s) = 2s 4 + 10s 
no 00 



b sx,± 



1 



2(s + 2)(s + 4) 
9s 2 - 22s - 5, 

(»-!)(« -2) 



(311) 
(312) 

(313) 
(314) 



(2s - l)(2s - 3) 
16si(s! + l)(s 2 -4s + l) + 6g(s) 



± 



2(s + 2)(s + 4) 



TO 

-7/o,i( s ) 

'/Jo,2( S ) 



-4s 5 - 19s 4 + 54s 2 - 2s - 5, 
(s-l)(2s-l) 
2s 

(s-l)(2s 2 -4s + 3) 
2s(2s- l)(2s-3) 



(315) 
(316) 

(317) 
(318) 



We note that a 

,1/2 



1/2,+ 
He, j 



(2) 



1/2, H 



(1) 



&r r ' (1) = for j = 1,2, as expected. However, by 
combining a/ro~o(s) anc ^ ^j/o ( s )' we nave 

1 16*i(ai + l)-d8(*) 



L IIa.2 



U IIo 



2 4(s + 2)(s + 4) 
d°(s) = 2s 3 + 13s 2 + 22s + 5, 



(319) 
(320) 



which differs substantially from the version with si 
1/2. In addition, it is elementary to combine cJ Io 2 (s) 



b°%-(s)}/2. We find 



-IIo,2 



U IIo 



J IIo 



(f - 1) 

2s 



(321) 



which differs from cJ Iol (s) by the factor 2s — 1. At 
s = 1, the only relevant odd s value for si = 1/2, 
these are equivalent. In addition, as for Type I and the 
even crossings of Type II, there is no simple relation be- 
tween the single-ion parameters aVr' ~o( s ) + 2bVj\ ' (s) and 



1 1 \ rm 



L v y{q) 



cr" = ±l 

^ Ss',s+a' (c m <T' s -(<T' + l)/2 
ff'=±l 

(^3,^3/^,2(9) 



ct"=±1 



M u y{q) 




(322) 



s' 13 ,s 13 +<t"/i/,i(<Z) 



/i7,o(g) = e (/+(«) + £i,» 13 , aM f- (<?) 



3in(g x a) sin^a) 



/£i(g) - 2(l-cos[a( g:c + g y )] 



x / ^Sl3 + (o"" + l)/2,S 2 4 

\ Si ,s 1 ,s 



/£ 2 (?) 



^/-(9K 2 



S + (<T' + l)/2,Si3,S 2 4' 

fvf (1) = 2 ( 1 - cos Kfe + g 3/ )]) 

/^S13 + ((T" + l)/2,S24 
^ LZ S 1 ,S 1 ,CT'cr"s + tT"(tT' + l)/2y ' 

/±(<?) = 1 + cos(q x a) cos(q y a) 

±cos(q z c)[cos(q x a) + cos(q y a)], 



(323) 

(324) 

(325) 
(326) 

(327) 
(328) 



7. Hartree INS functions 



The functions L v y{q) and M v y(q) in the self- 
consistent Hartree INS S g (q,u>) in the induction rep- 
resentation are given by 



where the A™, C™, 



are given by Eqs. 
respectively. 



r si,s 3 ,s j '~ r S 1 ,S 3 ,S 1 'lz,x, 



1(175) . (|23t) T (|238l) V and '^255 jl '-ip^ i l 
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